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ABSTRACT 

We analyse the relationships between galaxy morphology, colour, environment and stellar 
mass using data for over 10^ objects from Galaxy Zoo, the largest sample of visually classi- 
fied morphologies yet compiled. We conclusively show that colour and morphology fractions 
are very different functions of environment. Both colour and morphology are sensitive to 
stellar mass. However, at fixed stellar mass, while colour is also highly sensitive to environ- 
ment, morphology displays much weaker environmental trends. Only a small part of both 
the morphology-density and colour-density relations can be attributed to the variation in the 
stellar mass function with environment. 

Galaxies with high stellar masses are mostly red, in all environments and irrespective 
of their morphology. Low stellar-mass galaxies are mostly blue in low-density environments, 
but mostly red in high-density environments, again irrespective of their morphology. While 
galaxies with early-type morphology do always have higher red fractions, this is sub-dominant 
compared to the dependence of red fraction on stellar mass and environment. The colour- 
density relation is primarily driven by variations in colour fractions at fixed morphology, in 
particular the fraction of spiral galaxies that have red colours, and especially at low stellar 
masses. We demonstrate that our red spirals primarily include galaxies with true spiral mor- 
phology, and that they constitute an additional population to the SO galaxies considered by 
previous studies. We clearly show there is an environmental dependence for colour beyond 
that for morphology. The environmental transformation of galaxies from blue to red must 
occur on significantly shorter timescales than the transformation from spiral to early-type. 

We also present many of our results as functions of the distance to the nearest galaxy 
group. This confirms that the environmental trends we present are not specific to the manner 
in which environment is quantified, but nevertheless provides plain evidence for an environ- 
mental process at work in groups. However, the properties of group members show little 
dependence on the total mass of the group they inhabit, at least for group masses > 1Q^^M.q. 

Before using the Galaxy Zoo morphologies to produce the above results, we first quan- 
tify a luminosity-, size- and redshift-dependent classification bias that affects this dataset, and 
probably most other studies of galaxy population morphology. A correction for this bias is 
derived and applied to produce a sample of galaxies with reliable morphological type likeli- 
hoods, on which we base our analysis. 

Key words: galaxies: fundamental parameters — galaxies: structure — galaxies: statistics 
— galaxies: evolution — galaxies: clusters: general 

* This publication has been made possible by the participation of more 
than 100,000 volunteers in the Galaxy Zoo project. Their contributions are 

acknowledged at http://www.galaxyzoo.orgA^olunteers.aspx| . f E-mail: steven.bamford@nottingham.ac.uk 
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1 INTRODUCTION 

The discovery that galaxies may be naturally classified, simply by 
their visual appearance, into two principal types: spiral and ellip- 
tical, came even before these objects were firmly established as 
external to our own Galaxy tHubbleil922) . With further study it 
became clear that the differ ent galaxy types are very differently 
distributed throughout space i Hubble & Humasor]ll93lh . Over the 
intervening decades our understanding of galaxies has improved 
dramatically. Nevertheless, explaining the appearance of galaxies, 
and how this varies with position, remains a central concern. 

The visual appearance, or morphology, of a galaxy is an indi- 
cator of its current internal structure and kinematics, which in turn 
are a result of the galaxy's developmental history. Galaxies mainly 
comprise two structures, a spheroid and a disc. In the most basic 
terms, elliptical galaxies are simply a spheroid, whereas other mor- 
phologies generally comprise a central spheroid along with a larger 
disc, which often contains spiral arms. These disc and spheroid 
components appear to develop in very different ways. That the 
distribution of morphologies changes as a function of position in 
the universe suggests variability in these developmental processes. 
This may be the result of cosmological variations, or the interaction 
of galaxies with their surroundings. Galaxy morphologies therefore 
offer valuable information with which to construct and constrain 
theories of galaxy formation and evolution. 

The relationships between a galaxy's local environment and its 
colour and emission line strengths have been particularly well stud- 
ied recently. Both of these observables are related to star-formation 
history, and thus physical processes quite different from those that 
determine a galaxy's morphology. Comparing the behaviour of 
colour and morphology versus environment may thus provide pow- 
erful clues to the mechanisms through which galaxies have devel- 
oped into the population we see today. 

One reason for the current popularity of classifying galaxies 
by colour and, to a lesser extent, emission line strength, is the ease 
and accuracy with which these quantities can be measured in mod- 
em surveys, such as the Sloan Digital Sky Survey (SDSS). Tradi- 
tional morphological classification of galaxies, on the other hand, 
is extremely time consuming - requiring a human visual inspec- 
tion of each object. Studies have thus far been limited to sam- 
ples containing seve ral thousand galaxies iNakam ura et al.ll2004l : 
iFukugita et~Zf[2007h . although th e MOSES project (Mor phologi- 
cally Selected EUipticals in SDSS: ISchawinski et al.l2007h visually 
inspected nearly 50000 galaxies in order to identify a clean sample 
of ellipticals. 

Attempts at automatic morphological classifications have been 
made, with varying success. The most common, and arguably most 
useful, of these are simple measurements which quantify the radial 
light profile of a galaxy, such as concentration a nd Sersic index 
(e.g. iKauffmann etal]|2003l : iBlanton et al]|2003ij) . These quanti- 
ties measure the dominance of a spheroid over any disc component 
present. However, a spiral galaxy with a bright bulge or nucleus 
is still a spiral galaxy. Concentration is more strongly related to 
luminosity than morpholog y, with more lu minous galaxies having 
more concentrated profiles jCavazzi et alJi2000.) . These quantities 
are therefore not true measures of morphology in the traditional 
and most discriminating sense (as discussed in the next paragraph). 
Indeed, lyan der Wei (2008) finds very different mass and environ- 
ment dependences for galaxy type-fractions based on concentration 
versus those derived from an indicator more closely related to vi- 
sual morphology. 

There are several avenues to more sophisticated automatic 



morphologies, some which use physical insight or statistical 
methods to naturally classify galaxies, while others aim to di- 
rectly reproduce classifications by professional astronome rs (e.g. 
Lahav et al., 1995; Simard et al, 2002; Conselice 2003; Goto et al.l 



20031; iBall et al.l |2004 lAllen et al.l l200d Blakeslee etal] |2006| ; 
Conselicell2006l ; iBall et alj|2008l ; IPark et al.h2007 ). While a num- 



ber of these approaches are very promising, they do not yet provide 
a direct equivalent to traditional visual morphology, generally rely- 
ing instead on the correlations between true morphology and other 
parameters, such as luminosity, colour and concentration. The hu- 
man eye has consistently proved better than computational tech- 
niques at identifying faint spiral structure in noisy images, the ap- 
pearance of which is a primary indicator of morphology, along with 
the relative luminosity of bulge and disk. The presence and form 
of galaxy spiral arms has important physical implications. These 
patterns are related to density waves propagating around the disk, 
and indicate the dynamical state of the galaxy. They are also an 
important, though not exclusive, mechanism for inducing star for- 
mation, and thus provide information on the process by which a 
galaxy is currently forming it's stars. In addition, as the visibility 
of spiral structures may persist for some time after star formation 
ceases, their appearance gives an indication of the timescale of any 
decline in star formation. Finally, to have confidence in any auto- 
mated technique we must compare with a large number of visually 
classified objects that cover the full range of galaxy appearance. 

The Galaxy Zoo project was born out of the need for reliable, 
visual morphologies for a large sample of SDSS galaxies. For the 
reasons mentioned above, the presently available automated meth- 
ods were deemed insufficient for the task. Our approach was to 
enlist the public's help to visually morphol ogically classify all the 
galaxies in SDSS Data Release 6 (DR6, lAdelman-McCarthv et al.l 
,2008.) which were targeted for spectroscopy; nearly one million ob- 
jects. Further details of the Galaxy Zoo project, including its moti- 
v ation, desi gn and the initial stages of the data reduction, are given 
in lLintott e t al. ( 2008). 

A wide range of science is possible with the Galaxy Zoo 
dataset. One initial aim of the project, the investigation of a pop- 
ulation ofj^rej_bluej_star-^OT early-type galaxies, is consid- 
ered in lSchawinski et al] ( l2008h . Another early aim was to measure 
the statistical prop erties of spi r al gal axy 'spin' orientations, inves- 
tigated in detail in iLand et alj l20og) . This present paper makes a 
start in exploiting the Galaxy Zoo data to study the dependence of 
morphology on a host of galaxy properties. Here we concentrate on 
the dependence of morphology on local galaxy environment, how 
this compares with the dependence of colour on environment, and 
the role of stellar mass. 



1.1 Structure of this paper 

Before describing our own study, we set the scene with a brief dis- 
cussion of previous, related work in Sec. 11.21 Following this, in 
Sec. |2l we provide details of the data products and sample defini- 
tions J2.1| (. stellar masses J2.2b . visual morphologies J2.3I > and en- 
vironmental measures J2.4b that we employ in this paper. Though 
not essential for understanding our main results, an important com- 
ponent of this work is a quantification of the biases present in the 
morphological classifications. These are discussed whilst describ- 
ing our morphologies in Sec. 12. 31 but the full details are deferred to 
Appendix lAl 

The main results of this paper are presented in Sec.[3]& Sec.|4] 
Firstly, in Sec.jS] we consider in detail the local relationships be- 
tween galaxy visual morphology, environment and stellar mass. 
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Then, in Sec.|4l we compare these morphology relations with those 
derived using colour to divide the galaxy population, and explore 
the origins of the differences we find. Our conclusions are sum- 
marised in SeclH 



1.2 Previous studies of morpliology and colour versus 
environment 

The classic papers studying the relationships between visual mor- 
phology and environme nt in the local universe date from the 
1980's and early 199 0's (Dressier 1980b: Postman & Gelleilll984 
IWhitmore & Gilmor e 1991; Whitmore et al. 1993). The study of 
traditional, visual morphology for large samples of local galaxies 
has since stood still for the past decade, due to the great effort 
required to perform such measurements. In contrast, as discussed 
above, it is now straightforward to measure colours and simple 
structural parameters for large galaxy samples, and these quantities 
have therefore taken precedence. It is crucial that we link recent re- 
sults on large surveys, such as SDSS, back to the wealth of earlier, 
traditional morphology studies on which much of our understand- 
ing is founded. 

In order to efficiently sample a wide range of galaxy environ- 
ments, many studies have concentr ated on th e region s in and around 
rich galaxy clusters. In particular, lOressleil jl980ah measured the 
morphology of approxima tely 6000 galaxie s in 55 nearby, rich clus- 
ters. Studying this dataset. lPressleJ fl980b ') found a strong relation 
between local galaxy surface density and morphological type frac- 
tions. As local density increases, the spiral fraction was found to 
fall steadily, the SO fraction rises in a corresponding manner, and 
the elliptical fraction increases sharply at the highest densit i es. 

The same dataset was revisited bv lPostman & Gellej jl984) . 
who found that the morphology-density relation extended 
smoothly to galaxy group environments identified in the Cf A red- 
shift survey (Huchra et al. 1983) and by Huchra & Gellei| ( Il982h 
(see Fig. Both Dressier and Postman & Geller identified local 
density, rather than distance to the centre o f the nearest group, as 
being more closely related t o mor phology. IWhitmo re & Gilmore' 
( Il99ll) and I Whitmore etal] h993l) again reanalysed the Dressier 
(1980a) data, but contrastingly claimed that groupocentric distance 
was a better indicator of morphology than local galaxy density. 
These opposing results have yet to be resolved. We aim to address 
this with the Galaxy Zoo dataset in a forthcoming paper. 

Studying the dependence of an automated measure of mor- 
phology on both local de nsity and group ocentric distance in the 
SDSS Early Data Release. lOoto et al] ( 12003.) found both relations 
to be strong. Their most striking result is evidence that the fraction 
of galaxies with intermediate morphological types (early-type spi- 
rals, Sa-b) increases with galaxy density before falling at the high- 
est densities. Late-type spirals (Sc-d), on the other hand, steadily 
decrease in fraction with local density, whereas the elliptical frac- 
tion increases sharply at the highest densities. These findings were 
interpreted as a suggestion that multiple mechanisms were at work 
in shaping the morphology-density relation. 

Many mechanisms have been proposed for the transformation 
of spiral g alaxies to earlier m orphological types in dense environ- 
ments; see lBoselli & Gavazzil (|2()06l) for a thorough recent review 
of the proposed mechanisms and the evidence for them. The mech- 
anisms may be crudely divided into those which simply stop star 
formation, which indirectly affects morphology by removing the 
appearance of spiral arms and reducing the prominence of the disk, 
and those which affect the stellar kinematics of galaxies as well as 
leading to a cessation of star formation. While both categories of 



mechanism can form SOs from spirals, only the latter can create 
elliptical galaxies. 

With the launch of the Hubble Space Telescope, studies 
of morphology versus environ ment have focused on measur- 
ing the evolution with reds h ift dSmail et al.lll997l; lOressler et al.l 
TW; 'Poggianti et aL| 1 1 9991 ; iTreu et al.l |2003| ; ISmith et alj |2005|; 
Postman et al. 2005). Trends in morphological fractions versus en- 
vironment are found at z ~ 0.5-1 that are similar to those mea- 
sured locally. The principal difference is a reduced fraction of 
galaxies with SO morphology in distant intermediate- and high- 
density cluster environments, balanced by a higher spiral fraction. 
Neither elliptical galaxies nor the field population show significant 
evolution. The static nature of the massive elliptical populations 
has been confirmed by dedicated studies of luminous, red galaxies 
(LRGs), which find ver y little evolution i n their number density or 
properties since z ~ 1 dWake et al]|200 ^. 

As mentioned above, environmental trends in galaxy 
colours and star-f ormation rate s have received much at- 
tenti on recently dPimbblet et al 



Gom ez et alj 2003 ; Balogh et al 



2002 



2004 



Lewis et"ai] 
Baldrv et al 



2002 



2004 



2006 ; Kauffmann et al."2004VBl anton et aril2005l;IWeinmann et al.l 
2006; Blanton & Berlind 2007). The general result is that the 
fraction of blue (star-forming) galaxies decreases with local 
density, in favour of red (passive) objects. Interestingly, the 
properties of the individual galaxy sub-populations do not appear 
to change substantially, while their fractions vary greatly. Studies 
at intermediate redshift find that the fraction of blue, star-forming 
objects in clusters has decreased substantially since z ^ 0.5-1 
( Butcher & Oemleil [19841 ; iDressler & Gunn|[T993 ; IPoggianti et al.l 
19991). This is often associated with the decline of the cluster spiral 

population dCouch et aLlll998l) . 

The advent of the halo model dCoorav & Sh"etiill2002h has led 
to a focus on disentangling the properties of galaxies at the cen- 
tre of thei r dark matter halo and satellites orbiting within a larger 
halo (e.g. Ivan den Bosch et al. 2008). The models frequently as- 
sume that all satellite galaxies have ceased forming stars. How- 
ever, a recent halo-model analysis of the colour-marked correlation 
function has found that a non-negligible fraction of blue satellites, 
which depends on satellite mass, is required to reproduce SDSS ob- 
servations (Skibba 200^. Another interesting finding of such stud- 
ies is a confo rmity between the morphologies of sate llites and cen- 
tral galaxies dWeinmann et al.ll2006l ; lAnn et al.l2008t) . In this paper 
we do not attempt to separate central and satellite galaxies, but con- 
sider the population as a whole. 

Together, all these results suggest that more than one physical 
process is responsible for the observed dependence of morphology 
and colour on environment: one which only acts at early times and 
is responsible for the excess of ellipticals in dense regions, and on- 
going processes that prevent further star formation in ellipticals and 
which transform disc galaxies to early- type morphologies. How- 
ever, many of the details remain uncertain and speculative. Even 
so, our current models of galaxy formation are unable to simulta- 
neously reproduce many of the general features of the observational 
picture. 



2 DATA 

2.1 Basic galaxy properties 

The basis of the data used in this paper is SDSS DR6 
dAdelman-McCarthv et a"i]|2008h . For our main analysis we con- 
sider only those galaxies in the Main Galaxy Sample, extended 
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objects with rpctro < 17.77 dStrauss et alj 1200^ . In addition 
we only use galaxies with measured spectroscopic redshifts (with 
z > 0.01). Our sample is therefore incomplete in high density 
environments due to fibre collisions l lvon der Linden et al. 2 007,) ; 
the spec troscopic completeness in rich clusters is est imated to be 
- 65% ( ivon der Linden et al.ll20oilYoon et alj|2008h . Two SDSS 
spectrograph fibres cannot be placed closer than 55 arcsec. How- 
ever, the fibre assignments were based solely on target positions, 
with no consideration of other galaxy properties, and in cases where 
multiple targets could only have a single fibre assigned, the target 
selected to be observed was chosen randomly. Incompleteness due 
to fibre collisions is therefore independent of galaxy properties. As 
our analysis considers only trends in the fractions of objects versus 
environment, rather than absolute numbers, fibre collisions have no 
significant effect on our analysis. We refer to this sample as om full 
sample; it contains 565798 objects and is used in demonstrating 
the classification bias in Sec. l2.3l and deriving corrections for it in 
Appendix IaI 

We further restrict our main analysis to a redshift range 
0.03 < z < 0.085. We refer to this sample as our magnitude- 
limited sample, containing 192960 objects. The lower redshift limit 
ensures that the morphological classification bias correction is sta- 
ble, as explained in Sec. 12.31 and Appendix |A] The upper red- 
shift limit is a compromise between the number of galaxies in the 
sample, the luminosity range over which we are volume-limited 
{Mr < —20.17 mag), and the reliability of the local density esti- 
mates. 

In order to combine or compare galaxies across a range of 
redshifts we must account for the redshift dependent selection bi- 
ases. To remove selection bias from the analyses in this paper we 
restrict the galaxies considered to those that would meet our ap- 
parent magnitude, size and surface brightness criteria if they were 
located at the upper limit of the redshift range considered. As mea- 
sures of galaxy size and surface brightness we use the radius con- 
taining 50 per cent of the Petrosian flux, Rso, and the average sur- 
face brightness within this radius, p^o, all from the r-band imag- 
ing. Given the upper redshift-limit we adopt, z < 0.085, the Main 
Galaxy Sample limits (r < 17.77 mag, R^o 1 arcsec and 
1-150 ^ 23.0 mag arcsec^^), and our assumed cosmology (see be- 
low), we thus limit to the subsample of magnitude-limited sample 
galaxies with Afr < —20.17 mag, .Rso > 1.6 kpc and absolute sur- 
face brightness fi^o < —13.93 mag kpc~^. We refer to this sample, 
which contains 125923 objects, as our luminosity-limited sample. 
Figure IaTI (in Appendix lAt shows the redshift distribution for our 
luminosity-limited sample, compared with that of spectroscopically 
observed objects in the SDSS Main Galaxy Sample. This sample is 
used for much of the analysis in this paper, with the exception of 
when we wish to explore down to low stellar masses. In this case we 
use the magnitude-limited sample but either apply a bin-dependent 
upper redshift-limit to ensure we are complete in each mass-bin, 
or use a 1/Knax weighting. These mass-limited samples are fully 
defined in the following section. 

Our photomet ry is from the SDSS DR6 Ubercal 

dPadmanabhan et al] 1200 Sh . Where subscripts are omitted, we 
use model magnitudes for colours and Petrosian magnitudes 
otherwise. All magnitudes are on the AB zero-point system 
and corrected for Galactic extinctio n. Absolute magnitudes are 
determined using KCORRECT v4_l_4 jSlanton & Roweisll2007l) . 

Throughout this paper we assume a Friedmann-Lemaltre- 
Robertson- Walker cosmology with — 0.3, JIa = 0.7, Hq = 
70kms~^ Mpc~\ 



2.2 Stellar masses 

We determine stellar masse s from rpetro and (t t — r)modei using the 
relation given in figure 5 of lSaldrv et al.l ( |2006|) . which is calibrated 
on the spectrally-determined stellar masses of iKauffmann et al.l 
feoOS) and Glazebrook et al. (2004). Being based on a single 
colour, these stellar masses may be expected to be of limited ac- 
curacy. They display a 0.13 dex scatter around the spectrally- 
derived masses on which they are calibrated. However, this is com- 
para ble to the qu oted uncertainties on the spectral measurements 
(K auffmann et al. 2003) and the scatter a nd systematic differe nces 
between different spectral mass estimates jGallazzi et alj2005l) . On 
average, our stellar masses therefore have a comparable accuracy 
to spectrally-derived estimates. In any case, our colour-based es- 
timates provide a reliable way to rank galaxies, and thus examine 
trends, with respect to stellar mass. 

With the adopted limits of our luminosity -limited sample, 
z < 0.085 and Mr < —20.17 mag, we become incomplete for red 
galaxies below a stellar mass of ~ 1O^'' "^A^0. In order to explore 
a wider range when binning by stellar mass we must ensure a high 
level of completeness in each bin. We achieve this, only when con- 
sidering bins of stellar mass, by lifting the sample absolute magni- 
tude limit, i.e. using the magnitude-limited sample, but further lim- 
iting the redshift range of objects which contribute to each bin, such 
that < 0.1 per cent of galaxies in each stellar mass range are omit- 
ted due to their being fainter than the survey apparent magnitude 
limit. For example, in the bin with log(A1*/Al0) = 9.9-10.1, 
99.9 per cent of galaxies are bluer than {u — r) — 2.96, and hence 
have log( A4, /Lr ) < 0.37 in solar units (from the calibration of 
iBaldrv et alj2006l) . For log M* IMa^ = 9.9 this corresponds to an 
absolute magnitude limit of Mr < —19.18. The survey apparent 
magnitude limit of r < 17.77 and assumed cosmology thus require 
that we restrict objects in this stellar mass bin to 2 < 0.055. We 
also ensure completeness in terms of size and surface-brightness 
in the same manner as for the luminosity-limited sample, by limit- 
ing these quantities to the range of values observable at the highest 
redshift considered, z < 0.085. Given the low redshift limit of our 
magnitude-limited sample, z 0.03, we are able to study galax- 
ies with stellar masses down to log A^*/Al0 = 9.5 in statistically 
useful numbers. We refer to this as our binned mass-limited sample. 

Occasionally we will wish to produce a relation for a sample 
which is complete down to a low stellar mass limit, but for which 
limiting the redshift range as described above would leave us with 
too few high-mass galaxies. In this case we use the magnitude- 
limited sample, but weight each object by 1/Knax, where Vmax is 
the volume in which it would have been possible to observe that ob- 
ject given its absolute magnitude, size, and surface brightness, and 
given the survey limits on the corresponding apparent quantities. 
With this weighting we can limit to XogM,/ Mq > 9.8 without 
excessive uncertainty, and refer to this as our Vma.yi-weighted mass- 
limited sample. 

2.3 Galaxy Zoo morphologies 

The morphologies utilised in this paper are derived from classi- 
fications by over 80,000 members of the international public as 
pa rt of the Gala xy Zoo project. This project is described in detail 
bv iLintott et al.l (2008). Briefly, each galaxy received several, in- 
dependent morphological classifications, each by a different user. 
The four possible classifications were labelled as: 'elliptical', 'spi- 
ral', 'don't know', and 'merger'. The 'elliptical' class, in addition 
to containing galaxies with elliptical morphology, also contains the 
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majority of SO galaxies, as will be shown later in this section. We 
therefore refer to it henceforth as the early-type class. The merger 
class is mainly comprised of interacting pairs, generally with tidal 
features, which may or may not be expected to eventually merge. It 
also suffers from some degree of contamination from galaxy pairs 
that are overlapping in projection, but not physically related. This 
is of no consequence to the present paper, however, as only a small 
fraction (< 1 per cent) of objects are identified as a merger by a 
majority of classifiers, and these classifications are not specifically 
considered further herein. The spiral classification was subdivided 
into 'clockwise', 'anti-clockwise', and 'edge-on/unsure', referring 
to the direction and orientation of the spiral arms. This was primar- 
ily fo r use in studies of galaxy spins l lLand et al]|2008l ; ISlosar et aEI 
I2OO8I) . However, it also provides us with an indication of whether a 
galaxy was classified as spiral due to either noticeable spiral struc- 
ture or a disky, edge-on appearance. 

The median number of classifications per object is 34, with 98 
per cent of our full sample having at least 20 classifications each. 
These classifications were processed into raw 'likelihoods' for a 
galaxy being of a given morphological type, directly from the pro- 
portion of classifications for each type. We denote these as pd, p^p, 
Pdk and pmg for early-types, spirals, 'don't know' and mergers, re- 
spectively. 

A choice that must be made is how to go from the measured 
morphological-type 'likelihoods' to individual morphologies and 
type fractions. This choice depends upon the intended usage of the 
data. If one requires reliable morphologies for individual galaxies, 
or samples containing just one type, then assigning classifications 
to galaxies with type likelihoods above a certain threshold is appro- 
priate. The downside of this is that it results in a significant fraction 
of the sample being unclassified, these being objects with a spread 
of likelihood between two or more types. The choice of thresh- 
old is somewhat arbitrary, depending on the 'cleanness' and size 
of the required sample. The characteristics of samples defined us- 
ing Gala xy Zoo thresholded classifications are considered in more 
detail bv lLintott et"aLl ilOOA 

Alternatively, if one is interested in type fractions, then the 
likelihoods themselves can be used in the statistics. This retains 
more information than thresholding, and all galaxies can be in- 
cluded, but it cannot provide classifications for individual objects. 
Figure [T] shows the Galaxy Zoo type fractions as a function of 
redshift, determined by both thresholding (p > 0.8) and directly 
using the raw measured likelihoods. Assuming that there is negli- 
gible evolution in the galaxy population over the redshift interval 
considered, and that the survey contains a similar distribution of 
enviroimients at each redshift, then the true type fractions should 
be constant with redshift. As these two assumptions are expected 
to be reliable, any trends in the observed GZ type fractions with 
redshift may be attributed to a combination of two biases. The first 
is selection bias, due to variation in the size and luminosity distri- 
bution of galaxies in our magnitude-limited sample, resulting from 
the apparent selection limits. The second is classification bias, the 
tendency for otherwise identical galaxies viewed at different dis- 
tances to receive different morphological type classifications, due 
to signal-to-noise and resolution effects. 

Figure [T] also illustrates the difference in type fr actions when 
weigh ted or unweighted type likelihoods are used (see lLintott et al.l 
I2OO8I) . The effect of the weighting is to further increase the appar- 
ent fraction of early-type galaxies with increasing redshift. How- 
ever, the effect is small and only significant when thresholded clas- 
sifications are used. Throughout the rest of this paper we use the 
unweighted type likelihoods, without thresholding. 



As trends are clearly present in Fig.[T] the raw type likelihoods 
are, in general, biased. In Fig.|2]we attempt to remove selection ef- 
fects below given redshifts by limiting in absolute magnitude. The 
thin lines in this figure correspond to the raw type likelihoods. They 
still show strong trends with redshift, indicating that classification 
bias is present in the raw Galaxy Zoo data. 

For otherwise morphologically-identical galaxies, apparently 
fainter and smaller objects are more likely to be classified as early- 
type due to the diminished spatial resolution and signal-to-noise. 
However, this classification bias can be statistically corrected for, 
by applying an adjustment to the raw type likelihoods. These ad- 
justed likelihoods, Pci,adj andp^p^adj, may then be used to generate 
type fractions that are unbiased with respect to galaxy luminosity, 
size and redshift. The thresholding approach may also be based on 
these de-biased likelihoods, or limited to considering galaxies in 
regions of parameter space where the biases are shown to be small. 
The classification bias and our procedure for removing it are ex- 
plained in detail in Appendix|A] For the remainder of this paper the 
type fractions we consider are estimated by averaging the de-biased 
type likelihoods over bins of galaxy parameter space. 

The performance of our de-biasing procedure may be judged 
from the thick lines in Fig.|2] These indicate the type fractions after 
application of the de-biasing procedure. In the redshift range that 
is free from selection effects the de-biased type fractions are ap- 
proximately flat and, in particular, the ratio of early-types to spirals 
is constant. This indicates that the de-biasing procedure is working 
as anticipated. The downturn in early-type fraction at low redshift 
for the sample with the fainter limit is likely due to two effects. 
The small volume at low redshifts leads to a reduced occurrence 
of bright, and hence preferentially early-type, objects. In addition, 
the de-biasing procedure relies on the lowest redshift bins to estab- 
lish a baseline luminosity-size distribution against which higher 
redshifts are compared. Therefore, the lowest redshift bins cannot 
be entirely corrected by the procedure we adopt. To avoid this is- 
sue with our bias correction we limit our main analysis samples to 
z > 0.03. This conservative cut reduces the number of objects in 
our luminosity-limited sample by only 3 per cent. 

To give an impression of the how the galaxy images translate 
into morphological-type 'likelihoods'. Fig. [3] shows examples for 
Psp ~ (0.1, 0.5, 0.9) forfour different luminosity-size bins at fixed 
redshift. One can see that the appearance of spiral or disk features 
become clearer with increasing psp. 

Finally, we highlight that SO galaxies are not considered as a 
separate class in this work, due to the difficulty in discriminating 
them from the other classes and our initially conservative expecta- 
tions of the abilities of our classifiers. Face-on SOs are difficult to 
distinguish from elliptical galaxies, whilst edge-on SOs may appear 
the same as edge-on spirals. Even those who are highly experienced 
in morphological classification struggle with these distinctions. 

The Galaxy Zoo classification scheme specifies that objects 
with visible spiral arms or an edge-on appearance should be clas- 
sified as spirals, the remainder should generally be assigned to 
the early-type classification. The 'don't know' option is available 
for those who decide the image quality is insufficient for them 
to make the distinction, or for objects that do not appear to be 
galaxies. However, this resort was rarely chosen, partly because 
the Galaxy Zoo interface does not explicitly indicate the image 
resolution, and under-resolved galaxies of any type often appear 
to have elliptical morphology. Face-on or moderately inclined SOs 
will thus be found in the Galaxy Zoo early-type class. We can test 
this by examining the Galaxy Zoo type-likelihood distributions for 
objects classified by experts. For this we use the largest sample of 
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Figure 1. Morphological type fractions from out full sample of raw Galaxy Zoo likelihoods, plotted as a function of redshift. Classifications are based (left) 
on a likelihood threshold of p > 0.8, and (right) directly on the likelihoods themselves. The thick and thin lines corresponds to the weighted and unweighted 
likelihoods, respectively. "Unclassified" refers to all objects which do not meet the threshold of p > 0.8 for any class, whereas "don't know" refers to the 
likelihood of an object belonging to the "don't know" class, i.e. the fraction of times the classifiers clicked the "don't know" button. 
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Figure 2. Morphological type fractions, from our fiill sample of unweighted Galaxy Zoo likelihoods, as a function of redshift. The thick and thin lines 
correspond to the de-biased and raw likelihoods, respectively. Selection effects have been removed below (left) z = 0.088 and (right) z = 0.185 by imposing 
the magnitude limits given in each panel. These redshifts are indicated by the vertical dashed lines. The vertical dotted line in each panel indicates the low 
redshift limit (z ^ 0.03) applied to all our main analysis samples. Note that the ratio of de-biased early-type and spiral fractions is roughly constant over the 
redshift range for which selection effects are not present. 
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Figure 3. Example images illustrating the appearance of galaxies with 
different morphological-type likelihoods. Objects in the left, centre 
and right columns have psp ~ 0.1, 0.5 and 0.9, respectively. All 
objects are at a redshift z ~ 0.05, but each row corresponds to 
a different luminosity-size bin, from top to bottom: (MryRso) = 
(-22.625, 7.75), (-22.375, 5.75), (-21.875, 4.25), (-20.375, 2.25), 
corresponding to the white dots in Fig. IA6I These bins lie in the transition 
region between early-types and spirals in luminosity-size space, and are 
therefore not representative of our general sample, but provide a convenient 
comparison of otherwise similar objects with different morphological-type 
likelihoods. Each image is scaled to lOiJso on a side. 



morphologically-classified SDSS galaxies t hat attempts to distin- 
guish SOs from the other types, provided bv'Fukug ita et alj j2007l . 
hereafter F07). The early-type and spiral type-likelih ood distribu- 
tions a re shown in Fig.|5](this reproduces Fig. 8 from lLintott et al.l 
( l2008l) with the addition of histograms showing the de-biased type 
likelihoods used in this paper). We show distributions for the F07 
ellipticals, spirals, and SOs, as well as ambiguous objects that are 
classified as E/SO or SO/Sa. 

The majority of objects that F07 classify as E/SO or SO have 
high Pel and low psp (both raw and de-biased). The SO/Sa galaxies 
span a range of intermediate pei and Psp- We might expect edge- 
on SOs to be assigned to the spiral class in Galaxy Zoo. However, 
Fig.|4]shows that extremely few of the F07 E/SO, SO or SO/Sa galax- 
ies receive high psp . For the F07 sample, despite there being one SO 
for every four spirals, SOs contribute just 3 percent to the Galaxy 
Zoo spiral fraction. Either edge-on SOs identified by F07 appear 
more like ellipticals than spirals to the Galaxy Zoo classifiers, or 
F07 and Galaxy Zoo both classi fy edge-on S O s as s pirals, or per- 
haps a combination of the two. iLintott et alj j2008l) demonstrate 
that this is also true when thresholded classifications based on the 
raw type likelihoods are used. In any case, the uncertain distribution 
of true SOs between the classes henceforth referred to 'early-type' 



and 'spiral' should be borne in mind when considering the results 
of this paper. 

2.4 Measuring environment 

2.4.1 Local galaxy density 

We estimate local galaxy density by following exactly the detailed 
prescription in Baldrv et al. ( 2006), but extended from DR4 to 
DR6. The local density for a galaxy is given 
where div is the projected distance to the A'^th nearest neighbour 
that is more luminous than Air ~ —20 (with a small evolution 
correction, see below). In our analysis we use an estimate of lo- 
cal galaxy density, E, determined by averaging log Ejy for = 4 
and 5. In addition, the area ndj^ is modified close to a photomet- 
ric edge; and E is determined by averaging the density determined 
using spectroscopic neighbours (|c Az] < 1000 km s^^) with that 
determined using both spectroscopic and photometric neighbours 
(|Az;| < Zcrr.gs, the 95 per cent confidence i nterval, for galaxie s 
with only photometric redshifts determined bv lBaldrv et alj|200^ . 
The highly continuous imaging coverage in DR6 thus increases the 
number of galaxies with accurate local density estimates signifi- 
cantly beyond that expected simply from the increase in area since 
DR4. 

The density defining population is constructed independently 
of the other samples used in this paper. It comprises all galaxies 
with r < 18 in DR6 and with Air < Alr^m-nit — Q[z — zo), where 
A/r,iimit = — 20 and Q = 1.6, zp = 0.05 account for evolution as 
determined bv lBlanton eTZI j2003bl) . With these parameters, local 
density measurements are reliable for z < 0.085, the redshift to 
which the density defining sample is luminosity-limited. For den- 
sity defining galaxies without spectroscopic redshifts, Air is deter- 
mined using the redshift of the galaxy whose environment is being 
measured. We determine the potential range in log E that may re- 
sult from both edge effects and the use of galaxies with only pho- 
tometric redshifts. Galaxies for which the full range of potential 
log E is larger than 0.75 dex are rejected from the analysis sample. 
This amounts to ~ 10 per cent of objects, reducing our luminosity- 
limited sample to 113579 objects when we are considering local 
density. This sample has a log E uncertainty range of less than 0.4 
and 0.15 for 75 and 50 per cent of galaxies, respectively. For 25 per 
cent of our luminosity-limited sample the density is based on only 
spectroscopic galaxies far from survey edges, and thus the uncer- 
tainty range attrib utable to t hese issues is zero. 

As noted bv Ivan derW el et al. (2007) it may be more mean- 
ingful to limit the density-defining population in terms of stellar 
mass, rather than luminosity. However, as . van der Wei et al., show, 
while this may change the absolute values and distribution of local 
density, it does not significantly affect the ranking of galaxy envi- 
ronments, and thus would have no effect on our conclusions. 

2.4.2 Group properties 

In the standard ACDM cosmology, the properties of galaxies are 
thought to have a strong dependence on the characteristics of the 
dark matter halo in which they reside dCoorav & Shettj2002l) . Even 
without assuming this background model, the empirical clustering 
of galaxies into groups suggests a natural scale on which to mea- 
sure environment. The local density around a galaxy, as measured 
by the distance to its fifth-nearest neighbour, may fluctuate signifi- 
cantly and rapidly compared with the timescales for galaxy evolu- 
tion. Therefore, recent environmental influences on the properties 
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Figu re 4. Distribution of (to p) early-type and (bottom) spiral morphological type likelihoods for galaxies in our luminosity-limited sample with classifications 
from iFukugita et al] <2007l F07). Each panel contains objects classified by F07 as having the indicated morphological type. The dashed and solid lines 
correspond to the raw and de-biased type likelihoods, respectively. 



of a galaxy may not be reflected by the local density measured for 
that galaxy now. On the scale of galaxy groups, changes in environ- 
ment will generally be slower and smoother. It is unclear, however, 
which environmental scales play a significant role. Determining 
this is complicated by strong correlations between different envi- 
ronmental measures. 

To probe the influence of galaxy gro ups we utilise th e 
C4 catalogue of galaxy groups and clusters jMiller et al.ll2005b . 
The C4 catalogue was constructed by identifying galaxy over- 
densities simultaneously in both three-dimensional position and 
four-dimensional colour-space. This method greatly reduces the 
twin problems of projection effects and redshift-space distortions 
in identifying physical galaxy groups. Tests indicate that the C4 
catalogue is ~ 90 per cent complete and ~ 95 per cent pure above 
Mci = lQ^'^h~^ Mp) and within the redshift range considered 
here jMiller et al.ll20()a) . The method for estimating the virial ra- 
dius, and hence virial mass, for the C4 gro ups has been refi ned 
beyond those given in the original catalogue iTorki et alj2008h . 

At the time of analysis the C4 catalogue was only available 
for DR5. Whenever considering group quantities we thus limit our 
sample to objects with spectroscopy in DR5, though we continue 
to use DR6 measurements. We will employ two estimates of envi- 
ronment based on the C4 groups: the distance from a galaxy to its 
nearest C4 group, and the mass of the nearest C4 group. 

In order to determine the distance from a galaxy to its nearest 
group we require both the line-of-sight and projected distances be- 
tween them. A first approach is to convert redshifts, via an assumed 
cosmology, directly to line-of-sight distances between the galaxies 
and their nearest C4 group, dc4- However, this neglects the peculiar 
motions of galaxies. These motions are small for isolated galaxies, 
but increase in denser regions. They comprise two effe cts. F irst, 
there is a coherent infall in the vicinity of groups ( Kaiser. 19871) . but 
the velocities associated with this are generally low and produce a 
fairly small effect on the inferred galaxy distances. Secondly, there 
are the random motions of virialised group members. These virial 
motions are significant compared with cosmological velocities, par- 



ticularly at low redshift. The effect is such that galaxies which are 
truly group members are inferred to be at very different line-of- 
sight distances: the 'Fingers of God' effect. To correct for this one 
must identify galaxies that are likely to be group members, and put 
them at the same line-of-sight distance as their parent group. Our 
procedure for doing this is detailed in Appendix|B] resulting in cor- 
rected distances, Dci, which are additionally normalised by the 
virial radius of the nearest C4 group. 

When using groupocentric distances, we limit the galaxies 
considered to those which are closer to a C4 group than they are 
to an edge of the DR5 spectroscopic survey region. We also per- 
mit a galaxy's nearest C4 group to lie beyond the redshift limits 
of our magnitude-limited sample. We thus avoid the issue of clus- 
ters possibly located just beyond the survey and sample bound- 
aries. Recall that our redshift limits are 0.03 < z < 0.085, and 
so our magnitude-limited sample is restricted to redshifts where 
the C4 group cat alogue is reasonably complete (~ 90 per cent: 
iMiller et al.ll2005l) . 

Although they measure different quantities, the environmen- 
tal measures based on groups correlate strongly with local density. 
Figure |5] shows the relationship between local galaxy density and 
our two measures of distance from a C4 group. At large distances 
from a C4 group (> 10 Mpc or i?vir) galaxies are found in the full 
range of local densities, although predominantly these densities are 
low. However, approaching a C4 group the minimum local density 
steadily increases, such that all galaxies in the vicinity of C4 groups 
are also in regions of high local density. 

Correcting for the 'Fingers of God' effect clearly improves 
the correlation between distance to a C4 group and local density. 
There remain a large number of objects in regions of fairly high 
local density which are distant from a C4 group. There are sev- 
eral potential reasons for this: (i) even though density and group 
distance may be strongly related, the C4 catalogue forms an incom- 
plete census of galaxy groups; (ii) the 'Fingers of God' suppression 
is not fully effective, it has a reduced efficiency for objects at inter- 
mediate densities - for which we cannot distinguish between true 
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Figure 5. Local galaxy density versus (left) raw and (right) corrected distance to the nearest C4 group. Contours enclose 10, 25, 50, 75, 90 and \ 
galaxies in the luminosity-limited sample. 
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group members and nearby, but gravitationally unbound, objects - 
and is designed to be conservative and under-correct in such cir- 
cumstances, preferring incompleteness over contamination of the 
group members sample; (iii) the redshift interval over which lo- 
cal density is measured is large enough to allow nearby groups to 
influence the local density estimate for some isolated objects, but 
note that we do limit the local density uncertainty, see Sec. 12.4.11 
(iv) high local densities may occur in structures which do not meet 
the requirements to be selected as C4 groups, e.g. filaments, too 
few members, inhomogeneous member colours, etc. All of these 
issues are likely together responsible for the spread in local density 
at large distances from C4 groups and the bimodality visible along 
the line of (anti-)correlation between density and group distance. 

Close to C4 groups, within a few Mpc or virial radii, we expect 
a spread in local density due to substructure - azimuthal variations 
in local density. Even so, local density is confined to higher values 
than seen in the more typical environments. For substructures con- 
taining more than five density-defining galaxies the local density 
estimate will be well localised on the sky. However, low density 
measurements around groups may be precluded by the method's 
limited line-of-sight resolution. Nevertheless, the overall correla- 
tion between local density and C4 group distance must be a real 
effect. 

Local galaxy density and distance to a C4 group are thus use- 
ful alternative estimates of environment, with different caveats for 
their interpretation. The agreement of results when using either of 
these methods will therefore be a good sign that the results are ro- 
bust with respect to the general concept of 'environment' . Any con- 
trasts in the results from these two methods may point to different 
effects at work in classical groups and in more general dense re- 
gions. 

For galaxies within the influence of a C4 group, the mass of 
that group may well be expected to determine the magnitude of any 
environmental effect. The local density for members of C4 groups 
is shown in the left panel of Fig.[6]as a function of C4 group mass. 
As Fig.[5]suggests, galaxies within the virial radius of C4 groups are 



confined to high local galaxy densities. This is clearly seen by com- 
paring the histograms for cluster members and all galaxies along 
the side of the left panel in Fig. [6] However, there is almost no 
trend in the distribution of local galaxy density versus the mass of 
the group. 

The right panel of Fig. [6] shows the same plot as the left, but 
now containing all galaxies within a fixed physical radius, 3 Mpc, 
of a C4 group. A trend with group mass is now much more appar- 
ent. Comparing the left and right panels of Fig. |6] illustrates that 
the primary role of group mass is in determining the scale of the 
region, through the virial radius, over which local galaxy density is 
elevated, rather than the determining the absolute degree of density 
enhancement. 



3 MORPHOLOGY VERSUS ENVIRONMENT 
3.1 Local galaxy density 

As a first look at the Galaxy Zoo morphology-density relation, and 
a final demonstration of the removal of classification bias, in the 
left panel of Fig. [7] we plot the early-type fraction based directly 
on the raw type likelihoods, as a function of local galaxy density. 
This is shown for objects in four redshift slices. Note the offsets 
between the redshift slices: the average raw early-type fraction in- 
creases with increasing redshift, due to classification bias. How- 
ever, in each redshift bin a clear relation can be seen, with early- 
type galaxies favouring higher density environments. 

To produce a single, reliable, local morphology-density rela- 
tion we use the de-biased type likelihoods described in Sec. 12.31 
and derived in Appendix IaI The resulting morphology-density re- 
lations are shown in the right panel of Fig. |7] for the same four 
redshift slices as the left panel. Clearly the de-biasing procedure 
works well, as the relations from independent redshift ranges now 
match. We can thus combine galaxies over the full redshift range 
we are considering, 0.03 < z < 0.085. The morphology-density 
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Figure 6. (Left) points indicate local galaxy density versus C4 group mass for members of C4 groups with Dqa < IHvir- The solid and dashed Unes through 
the points indicate the median and quaililes of log S in bins of Mci- Along the left-hand side of the plot appear histograms comparing the log S distribution 
for C4 group members (solid) with that for all galaxies (dashed). The frequencies for the histogram of C4 group members have been multiplied by four to 
improve visibility. (Right) as the left panel, but now for objects within a constant 3 Mpc radius of each C4 group (Z)c4 < 3 Mpc). The additional dotted line 
indicates the median relation from the left panel. 




logio(S [Mpc-2]) logio(E [Mpc ^]) 

Figure 7. Early-type fraction versus local galaxy density for four redshift bins {z ^ 0.037, 0.051, 0.064, 0.078). The left and right panels show the fractions 
determined from the raw and de-biased morphological type likelihoods, respectively. 



relation for our full luminosity-limited sample, containing informa- 
tion on 113579 galaxies, is shown in Fig. [8] We determine the un- 
certainties on the mean type fraction in each bin from the standard 
deviation of the type likelihoods divided by the square-root of the 
number of objects in each bin. The reliability of these uncertainties 
has been confirmed by bootstrap resampling. In all plots we show 



2(7 uncertainty ranges. (The colour-density relation is also shown 
in Fig.[8]and subsequent figures, for which red galaxies are defined 
as those on the red side of the u — r versus stellar m ass bimodality 
using the divider determined bv lBaldrvetaLl ( l2006h . However, we 
defer discussion of this relation until Sec.|4l) 

Our morphology-density relation is extremely well defined. 
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Figure 8. Early-type (green, thick, solid line) and red (purple, thick, dotted line) fractions versus local galaxy density for galaxies in (left) our luminosity- 
limited sample and (right) our Vmax-weighted mass-limited sample, complete for stellar masses log(jVf*/A^0) > 9.8, and constructed by applying 1/Vmax 
weightings to our magnitude-limited sample (see Sec. |2.2) . The early-type fraction is determined from the de-biased type likelihoods, as described in Sec. |2.3| 
The shaded regions indicate the 2-sigma statistical uncertainties on each equally-spaced log E bin. Thin lines show the contribution to each relation attributable 
to variation in the stellar mass function with environment (see Sec. I3.1.U . The early-type fraction is discussed alone in Sec. |3] and compared with the red 
fraction in Sec.|4] 



It is consistent with a monotonic function, with a significant gra- 
dient over almost the whole range of galaxy densities that occur, 
0.1 < E < 25 galaxies Mpc~^. Even outside this range, there is 
no evidence for the relation flattening off completely. While more 
complicated representations could be acceptable, the data suggests 
that the relation is most appropriately described by a single, smooth 
function. 

To compare our morphology-density re lation to previous 
studie s, in Fig.|9]we reproduce the relations from lPostman & Gelled 
l ll984 hereafter PG84), whi ch contains a co mbination of data 
for galaxies in clusters, from iDressled (IT98o3) , and groups, from 
the Cf A redshift survey jHuchra etal.lll983h and lHuchra & Gelled 
( Il982l) . In PG84 the local density is measured in a rather different 
manner to that employed in this paper. However, to qualitatively 
compare our relation to these classic results we attempt to convert 
our measured local surface densities (S) to the local volume density 
scale ip) of PG84. We correct for the different depths of the den- 
sity defining populations (Mb < —17.5 versus AU < —20.0) by 
empirically determining the number density ratio given these selec- 
tion limits, as a function of local surface density. We convert from 
surface to volume densities assuming the overdensity to be spheri- 
cal. Finally we adjust to Ho ~ 100 km s~^ Mpc"'^. The resulting 
relation is shown by the thick green line in Fig.|9] 

As discussed in Sec. 12.31 SO galaxies are not explicitly sep- 
arated from the other types in the Galaxy Zoo data. However, the 
trend in the combined PG84 elliptical -i- SO fraction with local den- 
sity is very similar to our relation. A comparison of the type frac- 
tions suggests that the Galaxy Zoo early-type class contains the 
majority of the SO population, as also indicated by our earlier com- 
parison with F07 in Fig.|4] We defer the separation of the elliptical 
and SO populations to a future Galaxy Zoo project. 

The remaining differences between our early-type fraction 



versus local density relation and that of PG84 may be attributed to 
a number of issues, but the most important is the differing selection 
functions of the Galaxy Zoo and PG84 samples. As shall shortly 
be demonstrated, the morphology-density relation is strongly de- 
pendent upon the stellar mass (or luminosity) of the galaxies con- 
sidered, and thus different sample selections will naturally produce 
different relations. 

3.1.1 Stellar mass dependence 

A significant advantage of the Galaxy Zoo sample, in comparison 
with previous visual morphology catalogues, is its size. It is possi- 
ble to subdivide our sample in terms of various properties, and still 
retain sufficient galaxies in each subsample to make reliable infer- 
ences. In particular, we can consider morphology trends in narrow 
bins of stellar mass. 

It is now well known that the stellar mass function of galaxies 
varies with environment (Balogh et al. 2001; Zehavi et al. 200i). 
This can be, to some extent, theoretically understood by the biasing 
of the dark matter halo mass function varying with the large-scale 
density field jSheth & Tormenll 19991 : IMo & Whitdl2002l ; IMo et all 
l2004ldSkibba et al.l200 6^ although there are still many unsolved is- 
sues in connecting haloes with galaxies. In the left panel of Fig. 1101 
we show how the distribution of stellar mass depends on local 
galaxy density in our luminosity-limited sample. The distribution 
of stellar masses steadily shifts to higher masses with increasing 
local density. 

The mass function of galaxies is also expected to depend on 
morphologi cal type {M ailer 200 ^), as it varies wi th other galaxy 
properties l lCrotonet aL,2005i : IBaldrv et alj |2006|) . We shall see 
this clearly later in this section. In dense large-scale environments 
we expect a greater fraction of high-mass haloes, and such haloes 
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Figure 9. The thin hnes reproduce the morpho logical type fractions ver - 
sus local galaxy volume density from figure 1 o f lPostman & GelleiHl984h . 
From low to high, the groups of lines correspond to ellipticals (red), SO 
galaxies (orange), and ellipticals and SO galaxies combined (purple). In each 
group, the solid and dashed histograms represent data from CfA groups 
and clusters, respectively iHuchraetal. 1983), and the line indicates the 
converted cluster relation of Dressier 1 1980a). For comparison, the points 
joined by a thick green line show the Galaxy Zoo early-type fraction versus 
local su rface density relation fro m Fig. [8] converted to the volume density 
scale of lPostman & GelleJ il984 , as described in the text. 



are preferentially inhabited by galaxies with high stellar mass and 
early-type morphology. We must therefore consider whether the 
observed morphology-density relation may simply be a conse- 
quence of these effects. 

One way we can test this is by determining the morphology- 
density relation that would be present in our sample on the assump- 
tion that the early-type fraction is only a function of stellar mass. 
To do so, we directly measure the dependence of the early-type 
fraction on stellar mass for galaxies in our luminosity-limited sam- 
ple that reside in low-density environments (log E < 0). Then, 
for each bin of environment, we determine the morphological frac- 
tion predicted by applying this relation to the distribution of stellar 
masses in that particular environmental bin. The result is a relation 
between early-type fraction and local density that is due only to the 
empirical dependence of stellar-mass distribution on environment 
and morphology on stellar mass, without any direct dependence of 
morphology on environment. This is shown by the thin, green line 
in the left panel of Fig. [8] 

The variation of the stellar mass function with environment 
therefore accounts for ~ 40 per cent of the morphology-density 
relation in our luminosity-limited sample. However, this effect has 
been enhanced by the fact that this sample is limited in terms of 
absolute-magnitude. At low stellar masses all galaxies in the sam- 
ple are blue, and preferentially spiral. As we shall see later in this 
paper, such galaxies exist primarily at low densities. If, rather than 
an absolute-magnitude limited sample, we use a stellar-mass lim- 
ited sample, then the contribution to the morphology-density re- 
lation due to the varying stellar mass function is reduced. This is 
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Figure 11. Early-type (green, solid lines) and red (purple, dotted lines) frac- 
tions versus stellar mass for galaxies in low-density environments (log S < 
0). In order to fairly sample to low stellar masses we limit the redshift 
range contributing to each point to ensure completeness. The thick lines 
show the relation for our luminosity-limited sample, while the thin lines give 
the relation for our Vmax-^eighted mass-limited sample. Note that below 
log{Mt /Mq) = 10.3 the absolute magnitude limit of our luminosity- 
limited sample results in an absence of red galaxies, and hence a deficit of 
early-types. The shaded regions indicate the 2-sigma statistical uncertain- 
ties on each equally-spaced log A4 * bin. 



demonstrated using our Vmax-weig/jferf mass-limited sample (lim- 
ited to log(A1,/Al0) > 9.8 and with completeness ensured by a 
correction) in the right panel of [8] Figure [TT] shows the re- 
lationship of early-type fraction versus stellar mass in low-density 
environments for both our luminosity-limited sample and Vmax- 
weighted mass-limited sample. It is these relations, along with the 
dependence of stellar mass on local density as illustrated in the left 
panel of Fig. [TO] that are used above in determining the contribu- 
tion of the varying stellar-mass function to the morphology-density 
relation. 

For the V^ax-weighted mass-limited sample, only approxi- 
mately 32 per cent of the morphology-density results from the vari- 
ation of the stellar mass function with environment. The remainder 
must be due to changes in the early-type fraction with environment 
at fixed stellar mass. For both samples the morphology-density re- 
lation at low densities (E < 1 Mpc~^) is entirely explainable by 
the environmental variation in the stellar mass function. At greater 
densities it appears that some process must increase the early-type 
fraction at fixed stellar mass. 

To investigate this in more detail, we utilise our binned mass- 
limited sample to divide our sample into several complete bins of 
stellar mass (see Sec. |2.2t . We make these bins as narrow as reason- 
able given the uncertainties on the stellar masses. They thus have 
width 0.2 dex in stellar mass. We only plot alternate bins for clarity. 
The left panel of Fig.[T2]plots the morphology-density relation for 
several stellar-mass bins. We see that the relations are offset from 
one another, with the overall fraction of early-types increasing with 
stellar mass as expected from Fig[TT] As we have seen, some of the 
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Figure 10. The distribution of stellar mass in our luminosity-limited sample as a function of (left) local galaxy density and (right) groupocentric distance. The 
lines trace the 1, 5, 25, 50 (thick), 75, 95 and 99th-percentiles of the stellar mass distribution in bins of environment. The stellar mass distribution shifts steadily 
versus local density, while versus groupocentric distance most of the change in the stellar mass distribution occurs within 0.2R^\^. The dotted horizontal 
lines indicate a stellar mass of lO^" "^ A^0, below which our luminosity-limited sample becomes incomplete for red galaxies. 
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Figure 12. (Left) early-type and (right) red fraction versus local galaxy density for galaxies in selected narrow bins of stellar mass: log(A1 s,/ A4q) = 9.5-9.7, 
9.9-10.1, 10.3-10.5, 10.7-10.9, 1 1.1-1 1.3, in order from the lowest to highest average /^i, /red - The shaded regions indicate the 2cr statistical error on each 
log S bin. 



morphology-density relation is due to this effect, in combination 
with the varying stellar-mass function. 

However, a significant trend of morphological fraction with 
local density is present in each stellar mass bin. The morphology- 
density relation is not simply a product of a morphology-mass 
relationship and varying stellar-mass function; it exists even at 



fixed stellar mass. This confirms the recent result of Ivan der Wei 
(2008), whi ch is base d on an automatic measure of morphology 
(Blakeslee et al. 2006'), and highlights the difference between mor- 
phology and concentration; concentration shows little environmen- 
tal dependence at fixed stellar mass. 

Surprisingly, the change in early-type fraction with environ- 
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ment is very similar at each stellar mass, being approximately an 
increase of 0.2 in /d over the local density range probed (log E = 
— 1-1.5). For low stellar masses this corresponds to a doubling of 
the early-type fraction, while for high masses it is a fractional in- 
crease of around a third. 

There is a hint that for lower masses the morphology-density 
relation is flat in low-density environments. For the whole range 
of masses studied here, > 1O^'^A^0, a correlation between 
early-type fraction and local density is clearly present for E > 1 



Mpc 
sity. 



, but these is little evidence for a correlation below this den- 



3.2 Group properties 

In the previous subsection we considered the behaviour of mor- 
phological type fraction versus the density of neighbouring galax- 
ies. However, the relative velocities of galaxies in groups and clus- 
ters implies that local density may change rapidly, on timescales 
potentially shorter than those for morphological transformation. 
The average local density a galaxy experiences is related to the 
mass of the group halo it resides within. In addition, some pro- 
posed mechanisms for environmental galaxy evolution, such as 
gas starvation and tidal effects, relate directly to global halo prop- 
erties rather than the local galaxy density. Several studies have 
thus considered whether galaxy morphology is more fundamen- 
tally related to global environmental properties, such as distance 
to the centre of the nearest g alaxy group (groupoce n tric distance, 
Whitmore & Gilmord Il99ll: Iwhitmore et alj Il993l: lOomez et all 



20031 : iBlanton & Berlindll2007b . It would also improve our confi- 



dence in our results if the trends can be shown to be independent of 
the detailed way in which environment is characterised. 

To examine the influence of galaxy groups on morphological 
type fraction we use the C4 group catalogue. The range of local 
galaxy densities, E, covered in groups of varying mass, Mci, is 
illustrated by the left panel of Fig.[6] Clearly members of C4 groups 
reside in regions of high local galaxy density, but the distribution of 
E varies very little across the range of C4 group masses. In Fig.|5] 
we show the variation of E with distance from the centre of the 
nearest C4 group. 

The variation of morphological type fraction versus group 
mass for C4 group members is plotted in the left panel of Fig. 1131 
There is no evidence for a change in the early-type fraction over a 
group mass range of 10^'' - 10^^ A^©- Another estimate of group 
scale is the inte grated red luminosity of the group members (e.g., 
iLin et aLl l2006h . As a check we thus additionally show the type 
fractions as a function of the summed r-band luminosity of all 
group members with spectroscopic redshifts. There is a slight trend 
of increasing early-type fraction with increasing group luminosity. 
However, as versus group mass, this trend is small compared with 
the trends versus local density. Galaxies at all stellar masses con- 
sidered in this paper display similarly flat trends versus group mass, 
simply with offsets due to the dependence of morphology on stellar 
mass. 

In Fig. [14] we show the variation of morphological type frac- 
tion versus the distance to the nearest C4 group for our usual 
luminosity-limited sample and Vmax-weighted mass-limited sam- 
ple. A strong dependence is seen, with the fraction of early-type 
galaxies increasing for galaxies closer to a group centre. This 
dependence is particularly strong within the group viral radius, 
Dc4 ^ Rvir, although a slight trend appears to continue out to 
larger distances, particularly for the luminosity-limited sample. As 
with density, the relation appears consistent with a smooth function 



of groupocentric distance, rather than a broken line as advocated by 
lGotoetal.1 ( [20031) . 



3.2.1 Stellar mass dependence 

As was done in Sec. |3.1.n for local density, we can determine the 
contribution to the dependence of morphology on groupocentric 
distance from the environmental variation in the stellar mass func- 
tion. The stellar mass distribution versus groupocentric distance 
is illustrated in the right panel of Fig. [TO] There is a gradual in- 
crease in stellar mass with decreasing groupocentric distance, until 
Dc4 ^ 0.27?vir, within which the stellar mass distribution rapidly 
shifts to masses ~ 3 times larger. The exact radius at which this oc- 
curs may depend somewhat upon our 'Fingers of God' correction 
(see Appendix int. However, it is clear that galaxies in the cores of 
groups are typically more massive than those in the field. 

We determine the relation between early-type fraction and 
stellar mass at large groupocentric distances (log Dci > 1) in 
each sample, and combine this with the dependence of stellar mass 
on groupocentric distance. The result is a relation between early- 
type fraction and groupocentric distance that is due only to the 
empirical dependence of stellar-mass distribution on environment 
and morphology on stellar mass, without any direct dependence of 
morphology on environment. This is plotted as the thin green lines 
in Fig. [14] Any small trend at large groupocentric distances is at- 
tributable to the varying mass function. Below _Dc4 ~ 3-Rvir, how- 
ever, the early-type fraction increases faster than the mass function 
can explain. In the cores of groups, Dci ^ 0.27?vir, the quickly 
varying stellar mass distribution causes a similar rise in the early- 
type fraction, sufficient to account for most of the difference rela- 
tive to the field. However at intermediate groupocentric distances, 
the varying mass function is unable to account for any of the in- 
crease in early-type fraction. This clearly points to a separate pro- 
cess occurring throughout galaxy groups which increases the early- 
type fraction at a given stellar mass. 

We divide the morphology-groupocentric distance relation 
into bins of stellar mass in Fig. [15] Due to the smaller number of 
objects for which we have reliable groupocentric distances, the stel- 
lar mass range we can probe with good statistics is less than that for 
local density. Nevertheless, we can still consider over an order of 
magnitude in stellar mass. Each stellar mass bin displays a similar 
trend, with offsets due to the dependence of morphology on stel- 
lar mass. At fixed stellar mass there is a small, steady increase in 
early- type fraction at Dci ^ 3i?vir, but no evidence for a trend at 
larger radii. 

The trends in morphology with groupocentric distance are 
very similar to those with respect to local density. As explained 
in Sec. 12.41 these two ways of characterising environment are mea- 
sured by different methods and are subject to different interpreta- 
tional issues. The agreement in the relations of morphology versus 
both these measures of environment demonstrates that these trends 
are general and robust. 



4 COMPARING MORPHOLOGY AND COLOUR 

It is informative to compare the trends of morphology versus envi- 
ronment with those of colour versus environment. As discussed in 
Sec.[T] morphology is generally considered to be determined by a 
combination of a galaxy's dynamical state and star-formation his- 
tory. Colour is an indicator of a galaxy's recent (< 1 Gyr) star- 
formation history, with no direct dependence on the spatial dis- 
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Figure 13. Early-type (green, solid line) and red (puiple, dotted line) fractions versus (left panel) group mass and (right panel) integrated r-band group 
luminosity for members of C4 groups (Dc4, < l-Rvir) in our luminosity-limited sample. Only bins containing more than 50 galaxies from at least 5 groups 
are shown. The shaded region indicates the 2cr statistical error on each equally-spaced bin. The large circles with an'ows indicate the mean early-type (filled 
point) and red (open point) fractions for galaxies that are not members of a C4 group. The early-type fraction is discussed alone in Sec. [5] and compared with 
the red fraction in Sec.|4] 
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Figure 14. Early-type (green, thick, solid line) and red (purple, thick, dotted line) fractions versus distance to the nearest C4 group normalised by its virial 
radius, for (left) galaxies in our luminosity-limited sample and (right) our Vmax -weighted mass-limited sample, complete for stellar masses log {M * /M ) > 
9.8, and constructed by applying l/Vmax weightings to our magnitude-limited sample (see Sec. l2.2t . The shaded regions indicate the 2a statistical error on 
each equally-spaced log Dqa bin- Note that in this and subsequent plots the log Dqa axis is plotted with positive values on the left and negative values on the 
light for easier comparison with the plots versus S. The right-most log Dqa bin additionally includes all galaxies with lower log Z)c4- Thin lines show the 
contribution to each relation attributable to variation in the stellar mass function with environment (see Sec. 13. 2.0 . The early-type fraction is discussed alone 
in Sec.|3] and compared with the red fraction in Sec.|4] 
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Figure 15. (Left) early-type and (right) red fraction versus distance to tlie nearest C4 group normalised by its virial radius, for galaxies in narrow bins of 
stellar mass: \og{Mt / Mq) = 9.7-9.9, 10.1-10.3, 10.5—10.7, 10.9-11.1, in order from the lowest to highest average /d. The shaded regions indicate the 2cr 
statistical eiTor for each log Dc4 bin. For each stellar mass range, the light-most log Dc4 bin additionally includes all galaxies with lower log Dc4. 



tribution of its stars. Differences in the behaviour of colour and 
morphology with environment will thus give indications of the en- 
vironmental mechanisms at work. 

Several previous studies have addressed this issue using auto- 
matic morphology proxies. Those which consider structural mea- 
sures of morphology, such as concentration or Sersic index, find 
that these quantities are much less dependent upon environment 
than they are on galaxy mass (Kauffma nn et alj|2004l : iHogg et al] 
[2004). In contrast, colour and star-forming fr actions are strongl y 
dependent upon both environ ment and mass tealdrv et al.ll2006h . 
However, Ivan der Well \ 2008[) has shown that using a more so- 
phisticated morphology proxy, which accounts for the presence 
of 'bumpiness' in the surface-brightness distribution, leads one to 
conclude that morphology is more dependent on environment than 
mass. Different morphological and structural measurements thus 
lead to apparently opposing conclusions. With the Galaxy Zoo 
dataset we can address these issues using traditional visual mor- 
phologies, without recourse to automatic proxies. As discussed in 
Sec. [T] visual morphology is primarily concerned with azimuthal 
structure, i.e. spiral arms, which are not well characterised by prox- 
ies, and which contains different information to that provided by 
radial structural measurements. 

In order to compare how colour and morphology depend 
on environment we must adopt a measure of colour and a cri- 
terion for dividing objects into red and blue samples. We follow 
lBaldrvetalJ ( l2006l) . using the optimal, stellar-mass-dependent, di- 
vider they determine to separate galaxies based on the bimodality 
of the u — r colour distribution . This utilises the SDSS model mag- 
nitudes j Stoughton et al. II2OO2'), and therefore the colour is effec- 
tiv ely centrally w eighted. We simply class objects above and below 
the lBaldrv et alj divider as 'red' and 'blue', respectively. We do not 
attempt to work with likelihoods of whether individual galaxies lie 
on the red or blue sequence, tho ugh this is possible by assuming 
the Gaussian fits of lBaldrv et al.l 



Figure[8]plots both the early-type and red galaxy fractions ver- 
sus local galaxy density, for our luminosity-limited sample (com- 
plete for Mr < —20.17), on the left, and Vina.x-weighted mass- 
limited sample (complete for log{Mt / Mq) > 9.8), on the right. 
We see that the qualitative trends for morphology and colour are 
similar, but the morphology-density and colour-density relations 
are offset from one another. Averaging over our luminosity-limited 
sample the fraction of red galaxies which have early-type morphol- 
ogy is 67 per cent. However, the differences in the shape of the 
early-type and red fraction trends in Fig. [8] demonstrate that this 
fraction is a function of environment. It is clear that morphology 
and colour are not equivalent ways of classifying galaxies, and that 
they are not equally affected by environment. We shall explore this 
further below. 

Moving onto group properties, the left panel of Fig.|13|shows 
that, as was seen for early-type fraction alone in the previous sec- 
tion, there is no clear trend of red-fraction with group mass. The 
right panel indicates a small trend versus the summed r-band lu- 
minosity of group members, though this is small compared with 
the v ariation with enviro r mient and, as we shall see below, stellar 
mass. iBlanton & Berlindl ilOOj) find a strong variation of red frac- 
tion with group luminosity. However, this is driven by lower group 
masses than we consider here. At high group masses their re d frac- 
tions match ours. This also agrees with lPoggianti et al.l ( I2OO6) who 



find that for SDSS groups with > 10 A4q (assum ing the veloc- 
ity dispersion to mass conversion given by eqn. 2 of iBiviano et al.l 
2006) there is no mass-dependence of the star-forming fraction, 
while below this there is a clear dependence. Interestingly, they also 
show that the plateau of star-forming fraction in high-mass groups 
has only set in since z ~ 0.6. 

Earlier we showed that the early-type fraction is closely re- 
lated to the groupocentric distance. In Fig. [14] we show that this 
is also true for the red fraction. The fraction of red galaxies with 
early-type morphology varies strongly with groupocentric distance. 
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As for density, however, there is a considerable difference in the the 
red and early-type relations versus groupocentric distance, in terms 
of both an offset and the shape of the dependence. 

As was done earlier for morphology (see Sec. 13.1. it , we deter- 
mine the contribution to the relations between colour and environ- 
ment due to trends in the distribution of galaxy stellar masses versus 
environment. The relations between red fraction and stellar mass 
for low density environments in our data is shown in Fig. [TT] for 
both our luminosity-limited and Vina.x-weighted mass-limited sam- 
ples. The trends of colour versus environment expected from these 
relations, combined with the environmental dependence of stellar 
mass (Fig.llOt, are indicated by the thin, dotted lines in Fig|8l&[T4] 
for density and groupocentric distance, respectively. We find that, 
as for morphology, the varying stellar mass function is inadequate 
to explain the rapid increase in the red fraction at intermediate to 
high local densities, and at intermediate groupocentric distances. 
Again, this is plain evidence for an environmental process at work 
in groups. 

In low density environments we have seen that any 
morphology-environment relation can be explained by the depen- 
dence of the stellar mass function on environment. In contrast, 
for colour there is a strong indication of a residual trend versus 
both density and groupocentric distance. The slope of the colour- 
environment relations expected due to the varying stellar mass 
function is not as steep as that observed. This strongly suggests 
that, while morphology is unaltered, colour is significantly affected 
by environmental variation at low densities. Studies of galaxies in 
void environments have previously foun d colour, and star forma- 
tion, to vary with lo cal galaxy density l lRoiasetal.ll2004l , l2005l : 
ICeccarelli et al.l2008h . Figures[8l& ll4l clearlv demonstrate that this 
effect is in addition to that expected from the varying stellar mass 
function, and that it is not related to changes in the ratio of early- 
types to spirals. 

4.1 Stellar mass dependence 

The colour-density and colour-groupocentric distance relations for 
narrow bins of stellar m ass are shown in the right panels of Fig.ll2l 
& [15] respectively. As Baldrv et alj ( l2006h have shown, the red 
fraction depends strongly on environment at fixed stellar mass, par- 
ticularly for low masses (< IO^^A^q), for which the red fraction 
varies from < 0.2 to > 0.8 between low and high^density envi- 
ronments. This is in excellent agreement with Hain es et a i] 1 I2006I) 
who find that the relation between star-formation history and en- 
vironment is very different for giant and dwarf galaxies (roughly 
above and below ~ respectively). However, rather than 

a sharp division between giant and dwarf galaxies, we see a con- 
tinuous, though rapid, decrease in the sensitivity of red fraction to 
environment with increasing stellar mass. 

The dependence of red fraction on environment and stellar 
mass contrasts with the trends for early-type fraction shown in 
the left panels of these figures. The early-type and red fractions 
vary similarly with environment at high masses, but for low masses 
colour is much more sensitive to environment. The shapes of the 
environmental dependences also differ, with the morphology rela- 
tions displaying their largest gradients at the highest densities and 
smallest groupocentric distances, whereas the colour relations rise 
most rapidly at intermediate densities and groupocentric distances, 
and flatten off in the densest environments. 

As mentioned earlier for morphology alone, the agreement of 
the relations versus groupocentric distance and local galaxy density 
demonstrate the robustness of these trends in galaxy properties with 



respect to a general definition of environment. We will consider the 
detailed differences between the morphological and colour trends 
with respect to these two environmental measures in a subsequent 
paper. 

In this paper, we focus our investigation on the transforma- 
tion of galaxies due to environmental mechanisms. However, it is 
apparent that many massive galaxies are red and have early-type 
morphology independent of their environment. For some reason 
their star formation has ceased, and in many cases their stellar dy- 
namics have become dominated by random motions, without an 
obvious environmental cause. Some of these galaxies may be fos- 
sil groups (Jones et al. 2003), and hence only appear to be in low- 
density environments. However, in any case, such massive, isolated 
galaxies are likely to be centred in a gaseous halo, which should 
be supplying them with fuel for star formation. Preventing gas- 
cooling, and hence star formation, in massive galaxies remains a 
puzzle, although active galact ic nuclei (AGN) m ay play a sig nif- 
icant role ICroton et aLlliood : ISower et allbood ; ICroton & Farrail 
I2OO8I) . There is even a possibility that AGN activ ity could directly 
affect the morphology of galaxies jFan et al.l2008l) . However, AGN 
are only powerful enough to have a significant impact in massive 
galaxies. Furthermore, AGN activity is unlikely to depend directly 
upon environment, although in principle AGN activity may vary as 
a result of the same enviro nmental processes that affect star for- 
mation. iHaines et al.l i2007h investigate the relative importance of 
AGN versus environmental mechanisms for preventing star forma- 
tion in galaxies of different masses. 

It is often implicitly assumed that the red population corre- 
sponds to objects with early-type morphology, and that blue cor- 
responds to spiral. While this is correct for the majority of objects, 
making the assumption that these correspondences hold true in gen- 
eral can be highly misleading. The differences between the red and 
early-type fractions imply the existence of substantial populations 
of 'unconventional' galaxies: red spirals and blue early-types. Fur- 
thermore, the contrast between their dependences versus both envi- 
ronment and stellar mass require that these 'unconventional' popu- 
lations are not a constant contaminant, but vary strongly. One can- 
not, therefore, make inferences about the environment and mass 
dependences of morphology from a consideration of colour alone. 

We now examine the dependence of the red spiral and blue 
early-type populations on environment and stellar mass. 

4.2 Red spirals and blue early-types 

In the left panels of Fig.[T6]c&[T7]we plot the fraction of galaxies that 
meet our red spiral and blue early-type criteria as functions of local 
galaxy density and groupocentric distance, respectively. We con- 
tinue to derive morphological type fractions directly from the de- 
biased Galaxy Zoo type likelihoods. Here we simply modify these 
for 'red spirals' b y reducing the likel ihood to zero for those ob- 
jects bluer than the lSaldrv et al.l j2006h divider. Likewise, for 'blue 
early-types' we remove objects redder than the colour divider. 

The fraction of red spirals displays a clear peak versus both 
local density and groupocentric distance. In low density regions 
and far from C4 groups, red spirals constitute ~ 16 per cent of 
the local galaxy population brighter than Air < —20.17. This 
fraction almost doubles with increasing local density or decreas- 
ing groupocentric distance, rising to a peak of ~ 28 per cent at ~ 6 
galaxies Mpc^^ or ~ 0.4_Rvir. At higher densities or in the cores 
of groups the red spiral fraction declines sharply. 

The fraction of blue early-types simply diminishes steadily 
from ~ 12 per cent to ~ 2 per cent from the least dense to most 
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Figure 16. The number of red spirals (red, thick, sohd hne) and blue early-types (blue, dotted line) in our luminosity-limited sample versus local galaxy density 
(left) as fractions of our whole sample and (right) as fractions of all spirals and early-types, respectively. Also shown for comparison is the fraction of face-on 
red spirals (orange, thin, solid line), scaled to have the same mean fraction as the full red spirals sample, as indicated by the scale on the right of the plot. The 
shaded regions indicate the 2-sigma statistical uncertainties on each log S bin. 
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Figure 17. The number of red spirals (red, thick, solid line) and blue early-types (blue, dotted line) in our luminosity-limited sample versus distance to the 
nearest C4 group normalised by its virial radius (left) as fractions of our whole sample and (right) as fractions of all spirals and early-types, respectively. Also 
shown for comparison is the fraction of face-on red spirals (orange, thin, solid line), scaled to have the same mean fraction as the full red spirals sample, as 
indicated by the scale on the right of the plot. The shaded regions indicate the 2-sigma statistical uncertainties on each log Dq^ bin. The right-most log Dq^ 
bin additionally includes all galaxies with lower log Dq^. 



dense environment s. The extreme members of this population are 
studied in detail bv lSchawinski et al] ( l2008h . and are also found to 
prefer entially inhabit low-density environments. ISchawinski et al.l 
( l2007h have also previously shown that the fraction of blue early- 
types declines with local galaxy density. The blue colours of these 
galaxies indicates that they have recently formed significant num- 
bers of new stars, in contrast to the usual pict ure of early-types 
as 'red and dead'. By using UV— optical colours. ISchawinski et al.l 
( l2007h have shown that over 30 per cent of bright (Mr < -21.5) 



early-types have recently formed stars. Our fractions are lower, due 
to the lower sensitivity of the (u — r) colour used in this paper to 
low levels of recent star formation, but still substantial, and trace 
the same environmental dependence. 

For both red spirals and blue early-types there is little evidence 
for any dependence on groupocentric distance beyond ~ 5i?vir- As 
previously, we find very little dependence of group members on 
the mass of the group in which they reside, over the range 10^^ '* - 
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The fraction C4 group members which are red spirals 
or blue early-types is flat at ^ 26 and ~ 6 per cent, respectively. 

What we have discussed so far, and plotted in the left panels 
of Fig. [16] & [T71 is the fraction of all galaxies that are either red 
and spiral or blue and early-type. This quantity indicates the envi- 
ronments in which these objects contribute most significantly to the 
overall galaxy population. However, it is also informative to con- 
sider the fractions of spirals that are red, and early-types that are 
blue, as the spiral and early-type fractions are themselves strong 
functions of environment. These 'type normalised' fractions are 
plotted in the right panels of Figs. [T6] & [TT] as functions of local 
density and groupocentric distance, respectively. 

The environmental dependence of the fraction of spirals that 
are red is simpler than the fraction of the whole galaxy population 
that are red and spiral. We now see that the decline in the fraction 
of red spirals at high densities, and in the the cores of groups, is a 
result of the declining spiral fraction. The fraction of spirals that are 
red increases steadily to the highest local densities. The dependence 
versus groupocentric distance is most interesting, as it appears lo- 
calised to a range of radii. Between 0.4 and 6i?vir the fraction of 
spirals that are red doubles, while outside of this range it is con- 
stant. A similar behaviour is seen for the fraction of early-types 
that are blue, with a two-thirds decline over the same groupocen- 
tric distance range. 

To gain further insight into the behaviour of the red spiral and 
blue early-type populations with environment, we again utilise our 
binned mass-limited sample and split our sample into complete, 
narrow bins of stellar mass. Fi gures [Tsl & 1 1 9 1 plot the fractions of 
spirals that are red and early-types that are blue for these stellar 
mass bins, versus local galaxy density and groupocentric distance, 
respectively. 

At all stellar masses there is a trend in the fractions of spi- 
rals that are red and early-types that are blue versus environment. 
In low density environments these fractions are strong functions of 
stellar mass. In these environments, at high stellar masses, A^, > 
10^^ Mq, the majority of spirals, and almost all early-types, are 
red; at low stellar masses, Mt, < lO^'^A^o, almost all galaxies 
are blue, irre spective of their morphology. This corresponds with 
the finding of' Conselicd ( l2006h that the correlation between galaxy 
colour and stellar mass does not strongly depend on morphology. 
It does, however, depend strongly on environment. The environ- 
mental dependence for low mass galaxies is much stronger than for 
high mass galaxies, such that in high density environments (E > 10 
Mpc~^, Dc4 < O.Siivir) 80 per cent of spirals and 90 per cent of 
early-types are red, irrespective of their stellar mass. In dense en- 
vironments the fractions of spirals that are red and early-types that 
are blue are therefore only weak functions of stellar mass. 



4.2. 1 SO galaxies and red spirals 

We see a significant fraction of red spirals in the field, particularly 
at high stellar masses. This results from a lack of correspondence 
betwe en the morphology and colour bimodalities jvan den Berghl 
l2007h . which is present across all environments. Many massive 
galaxies appear on the red side of the colour bimodali ty, irrespec- 
tive of their morphology. As shown in appendix A of Ijames et al.l 
( l2008h . with typica l selections one fin ds that half to two-thirds of 
Sa galaxies are red. lMaller et al. iiooi) have also shown that many 
red (and concentrated) objects are truly disk galaxies. These would 
often be mistakenly considered to be early-types in studies using 
colour and concentration as proxies for morphology. About one- 



third of these objects are red simply due to being orientated edge-on 
and thus suffering from significant dust reddening. 

A substantial fraction of the red spiral population we have 
identified shows a strong dependence on environment, preferen- 
tially occurring in intermediate local densities and within the virial 
radius, but not the cores, of galaxy groups. It is unclear whether 
our red spiral sample comprises only similar objects, which dis- 
play only limited environmental variation, or whether it is made up 
from several distinct populations with different environmental de- 
pendences. However, we do not see any change in the fraction of 
red-spirals that are classified as edge-on/unclear as a function of 
environment. Thus, the environmental trends cannot simply be at- 
tributed to a change in the proportion of red-spirals that are selected 
without having clearly visible spiral arms. 

SO galaxies possess discs and are generally found to have lit- 
tle star formation and to be red in colour. They are also known 
to preferentially inhabit clusters, especially at low redshifts (e.g., 
lDressleilll980bl) . This morphological type is not explicitly identi- 
fied in the current Galaxy Zoo classification scheme, and SOs are 
thus distributed between the early-type and spiral classes in an un- 
certain manner. There is therefore a concern that the environmental 
behaviour we have found for red spirals may simply be due to the 
SO population. However, recall that Fig.|4]demonstrates that almost 
all the objects classified as E/SO or SO by F07 have low spiral likeli- 
hoods in Galaxy Zoo. These objects therefore contribute very little 
to the spiral fraction. Objects classified by F07 as SO/Sa also gen- 
erally have low spiral likelihoods, and many of these, especially 
those with higher spiral likelihoods, may be better classified as spi- 
ral rather than SO, as they must have a suggestion of spiral arms. It 
is thus apparent that the red spirals we discuss above are an addi- 
tional population to the SOs considered in many other studies. 

There remains the possibility that Galaxy Zoo, F07, and thus 
presumably most other studies, mistakenly classify a significant 
fraction of edge-on SOs as spirals. In Galaxy Zoo we can ask the 
question of whether objects were classified as spiral due to the pres- 
ence of visible spiral arms or simply due to a disky, edge-on appear- 
ance. We define as 'edge-on/unclear' those objects for which the 
majority of Galaxy Zoo classifiers could not individually discern 
a spiral arm direction {pew + PACw < Peu where CW, ACW 
and EU refer to clockwise, anti-clockwise and edge-on/unclear, re- 
spectively). Note that this is conservative, as if even a minority of 
classifiers can discern spiral arms in a galaxy then it is very likely to 
be a spiral. The classification bias correction (see Appendix [At in- 
creases the spiral likelihood for such objects when they have appar- 
ent luminosities or sizes which may make distinguishing them from 
early-types difficult. The spin likelihoods (pcw , Pacw ,Peu) have 
not been corrected for such biases, and simply reflect the propor- 
tion of user classifications. Spiral galaxies may therefore include 
a surprisingly high fraction of edge-on/unclear objects, given the 
above definition. 

For all objects that are more likely to be spirals than anything 
else (Psp.adj > 0.5), the fraction of edge-on/unclear objects is 55 
per cent, while for red spirals it is 81 per cent. There are three pos- 
sible reasons for this discrepancy: (a) edge-on spirals suffer from 
greater dust-reddening, so will be more prevalent in a red sample; 
(b) due to the connection between star formation and spiral arms, 
objects with low specific star-formation rates will tend to be red 
and have lower spiral arm contrast, and are thus more likely to be 
classified as edge-on/unclear; and (c) a fraction of objects classified 
as edge-on/unclear spirals are actually SO galaxies, with typically 
redder colours. For cases (b) and (c), there must be a population 
of these objects viewed closer to face-on and hence classified as 
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log^o(S [Mpc ']) logio(S [Mpc ^]) 

Figure 18. The fraction of (left) spirals that are classified as red galaxies, and (right) early-types classified as blue galaxies, plotted as a function of local galaxy 
density, for our binned mass-limited sample: log(A4,/A^Q) = 9.5-9.7, 9.9-10.1, 10.3-10.5, 10.7-10.9, 11.1-11.3. The shaded regions indicate the 2-sigma 
statistical uncertainties on each log S bin. 
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logio(Dc4 [Rvir]) loglo(Dc4 [^vir]) 

Figure 19. The fraction of (left) spirals that are classified as red galaxies, and (right) early-types classified as blue galaxies, plotted as a function of groupocen- 
tric distance, in our binned mass-limited sample: log{A1*/A4Q) = 9.7-9.9, 10.1-10.3, 10.5-10.7, 10.9-11.1. The shaded regions indicate the 2-sigma 
statistical uncertainties on each log Dq^ bin. 



early-type. The environmental trends we see in a fraction of the 
red spirals must correspond to some combination of environmental 
variations in the (a) dust-content, (b) star-formation rate, and (c) SO 
contamination of the spiral galaxy population. 

Edge-on/unclear spirals are, in the median, 0.29 mag redder 
than galaxies with clearly visible spiral arms (i n (u — r)niodei. rel - 
ative to the stellar-mass-dependent divider of iBaldrv et alj|200^ . 



However, determining whether this is due to extinction affect- 
ing otherwise normal spiral galaxies, or due to older luminosity- 
weighted stellar populations, is beyond the scope of this paper. 

For 19 per cent of our red spirals significant spiral arms are 
seen by a majority of classifiers, and the true fraction with spiral 
arms, including those undetected in Galaxy Zoo, must be higher. At 
least a fifth of our red-spiral population must therefore be true spi- 
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rals, not SO galaxies. We can check that these face-on spirals, which 
we are certain possess spiral arms, behave in the same manner as 
our whole red spiral population. This is shown by the orange lines 
in Figs.[T6]&[T7] which plot the fractions for face-on red spirals 
only, scaled by the ratio of all red spirals to face-on red spirals. As a 
function of local density, the fraction of face-on red spirals behaves 
very similarly to the full red spiral sample. Versus groupocentric 
distance they also behave similarly, except within Dqa < 0.47?vir, 
where there is an excess of face-on red spirals compared with the 
full red spiral population. This suggests that in groups edge-on red 
spirals are less likely to be classified as spirals than they are in the 
field. This is interesting in itself, but here we simply stress that this 
implies that our red spiral population is not dominated by SOs in 
groups, rather there may be proportionally fewer SOs in our red spi- 
ral sample in dense regions. A substantial fraction of our red spiral 
population thus appears to comprise galaxies with significant spiral 
arm structure, even if many are seen edge-on. 

As explained in Sec. 12.11 the colours used in this study are 
based on SDSS model magnitudes. They are thus sensitive to the 
dominance of the bulge component of a disc galaxy. In contrast, 
the Galaxy Zoo morphologies are sensitive to the presence of spi- 
ral arms or an edge-on disk, even when the galaxy light is dom- 
inated by a bulge. This is because the images were displayed for 
classification using a non-linear (arcsinh) stretch, as described in 
iLupton et al] ^004). Other photometric apertures, such as Pet- 
rosian, may give colours that are more representative of the over- 
all galaxy, and so which may be more suitable for identifying 
galaxies with red discs. Measuring colour in an annular ap erture 
may be most appropriate for this purpose jPark et ^l2007t) . The 
choice of other photometric bands (e.g., g — r) and dividing lines 
in the colour-magnitude diagram may somewhat modify the trends 
shown. However, we do not expect the overall behaviour of the red- 
spiral population versus environment to depend dramatically upon 
our choice of colour. 



4.2.2 Breaking down the dependence of red fraction on 
environment 

Having studied the relationship between morphology and environ- 
ment, and the environmental variations in the fractions of early- 
types and spirals that are red, we are now in a position where we 
can identify the various contributions to the dependence of red frac- 
tion on environment. 

The red fraction can be written as a sum of contributions from 
red spirals and red early-types: 



r /red — sp r i / 1 ./blue — cl i <. 
/red — J ,/sp + i J I /ol ■ 



/sp 



.fol 



(1) 



We have already determined /sp, /d, /rod-sp//sp and 
/biuc-Gi//ci, as functions of density. To assess the red frac- 
tion resulting from the morphology-density relation we simply set 
/rcd-sp//sp and /biuc-ci//ci to their average values in Eqn. [T] 
while allowing /sp and /d to vary with density. The contribution 
from the changing fraction of spirals that are red is found by 
holding all terms except /rcd-sp//sp constant at their mean values. 
Finally, by holding all terms except /biuG-ci//oi at their mean 
values we determine the contribution from the environmental 
variation in the fraction of early-types which are blue. 

These contributions to the relation between red fraction and 
local density, for our luminosity-limited sample, are plotted in the 
left panel of Fig. |20l In order to examine only trends, rather than 



absolute values, we normalise each curve to the mean value of the 
three lowest density bins. The most striking result of this exercise, 
though already partly clear from Fig. [8] is that the morphology- 
density relation is responsible for very little of the colour-density 
relation. The modest increase in the early-type fraction with lo- 
cal density is not enough to explain the much larger increase in 
red fraction. Furthermore, as the early-type and spiral populations 
both contain significant fractions of red galaxies, which increase 
to comparable proportions in dense environments, the rise in the 
early-type fraction with density has surprisingly little effect on the 
red fraction. The colour-density relation is driven by changes in 
the red fraction at fixed morphology, i.e. the fraction of early-types 
and spirals that have red colours. In low density environments the 
rise of red spirals and fall of blue early-types contribute similarly. 
At higher densities the rapidly increasing red spiral fraction is the 
dominant contribution to the colour-density relation. 

Note that the variation of the stellar mass function with en- 
vironment is still present in Fig. |20] Its influence will be shared 
amongst the various contributions, for as we have seen, early-type 
fraction and the fraction of early-types and spirals that are red de- 
pend strongly on stellar mass (Fig.[T2]and Fig.llSt. 

The relative importance of the red spiral population varies 
with stellar mass. To show this we determine the contribution to the 
colour-density relation over the range log E = 0.8-1.5 in narrow 
bins of stellar mass. These are shown in the right panel of Fig. 1201 
Using narrow stellar mass bins also has the advantage of remov- 
ing the effect of environmental variations in the stellar mass func- 
tion. Red spirals dominate the colour-density relation at low stellar 
masses, but at high masses red spirals, blue early-types, and the 
morphology-density relation contribute equally. 

4.2.3 Why are galaxies red in dense environments? 

Our red spirals, and particularly the portion which are prevalent in 
the outskirts of galaxy groups, bear a strong resemblance to popula- 
tions discussed by a number of previous studies. These include the 
anaemic spirals of van den Berghl l ll976l) . the passive spira ls iden- 
tified by iGoto et afl ( l2003h and confirmed by [V amauchi & Goto! 
('2004'), Ha truncated (e.g.. iKoopmann & Kenne v 2004) and HI- 
deficie nt spirals (e.g.. IVogt et al.l 2004h . More recentlv Lane et al.l 
( l2007h have studied a population of dusty, red spirals associated 
with the i ntermediate density regions in a supercluster, which 
IWolf etalj ( 2007) and Wolf et al. ( 2008) find to still be hosting low 
levels of highly-obscured star formation. All of these populations 
are obviously closely related, and their memberships must over- 
lap considerably. However, the various selections may highlight 
objects at different stages of an evolutionary sequence, or subject 
to the action of different mechanisms. In general, most of these red 
spiral populations appear to result from the removal of gas from the 
haloes and/or discs of normal spiral galaxies in cluster and group 
environments. 

Considering Fig. [16] one might conclude that a more general, 
rather than a group-specific, mechanism is responsible for the envi- 
ronmental variation in the red spiral fraction. This is implied by 
the finding that the slope of the red-spiral fraction versus local 
density relation declines smoothly and does not entirely flatten- 
off at low densi t ies. Fo r example, such a conclusion is reached 
by iBalogh et Zl (l200il due to a similar behaviour of the star- 
forming fraction versus local density. However, when considering 
groupocentric distance as the environmental indicator, in Fig. [17] 
we see that almost all the variation in the red-spiral fraction occurs 
within a few i?vir. The mass of the group involved does not signif- 
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Figure 20. (Left) the change in the red fraction of galaxies in our luminosity-limited sample, divided into contributions due to the variation with local density 
of the morphological type fractions (green, dotted line), the fraction of spirals that are red (red, dot-dashed line) and the fraction of early-types that are blue 
(blue, dashed line). The total colour-density trend is shown by the solid, purple line. (Right) the contiibutions to the change in red fraction over the local 
density range log S = 0.8-1.5 as a function of stellar mass (indicated by the vertical dotted lines in the left panel). 



icantly affect the red-spiral fraction, at least above some threshold 
which lies below the range of group masses considered in this paper 
(> lO^TWo). 

Recall from Fig. |6] that there is a spread of local density for 
group members: and thus a sharp transition in groupocentric dis- 
tance will result in a smoother variation with respect to local den- 
sity. This argument suggests that the majority of the environmental 
dependence of red-spirals is due to group membership, rather than 
local density, and hence that a significant proportion of red spirals 
are a result of group-specific mechanisms. 

In the cores of groups, many galaxies have early-type mor- 
phology and red colours. Much of this can be explained by a 
shift of the stellar mass function to higher masses in these re- 
gions, together with the higher early-type and red fractions for 
more massive galaxies (Figs [TO] & lilt . However, this does not 
help to explain the high early-type and red fractions in the out- 
skirts of groups. At a given stellar mass the early-type fraction 
is higher in dense environments. This indicates that some pro- 
cess must have created early-types preferentially in dense regions, 
either at the epoch of galaxy formation or subsequently. A re- 
ce nt direct study of evoluti on in the morphology-density relation 
by Ivan der Wei et al.l ( l2007h demonstrates that the the relation has 
not significantly changed out to z ~ 1, when one considers a 
mass-selected sample. The evolution f ound by previous studies 
dSmith et aLll2005l : 1Postman et al.ll2005t) is attributed to effects of 
their luminosity-based selection. The star-formation properties of 
galaxies in dense regions appears to have changed, while their 
morphologies have not. This i s at odds w ith the picture of spirals 
transform ing to SOs in clusters jPressler et al. 1997; Bedregal et al 



ative to those of 
7! '" 



spirals (iBursteiij [19791: iDressleJ 1 1 980bl Ig 
: Zabludoff 2004 : BoseUi & GavazzilllOm 



2006; Ar agon-Salamanca et ai]|2006l ; Isarr et al.ll2007l ; lDesai et al 



2007) . However, this picture remains contentious, with other stud- 
ies concluding that luminous SOs cannot have formed from faded 
spirals because of their larger bulge sizes and luminosities, rel- 
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1980; Christlein & Zabludoff 2004"; BoseUi & Gavazzill2006h . The 
study of van der Wei et al. (2007) is limited to massive galaxies 
{M, < IQ^'^-^Mq), so there is still room for the morphology- 
density relation to evolve at lower masses. 

The significant population of red spirals we identify, which 
preferentially reside in the outskirts of groups, are predominantly 
of intermediate stellar mass. Furthermore, the lowest mass objects 
appear to be those which are most strongly affected by their envi- 
ronment, as in the field they are almost all blue. The star formation 
in these objects must have declined recently, within the past few 
Gyr, in order for their spir al arms to r emain with sufficient con- 
trast to be visible. Indeed, IWolf etalj ( l2007h find that when one 
considers stellar population age, rather than simply the colour bi- 
modality, the morphology-density relation may be explained by 
an age-density relation and age-morphology relation (using more 
finely divided morphologies than our study). Their results imply 
that, in clusters, red spirals are younger than red early-types, and 
that both the morphology of a galaxy and its distance to the cen- 
tre of the cluster are strongly related to the length of time since 
the galaxy's star formation ceased, presumably when it entered the 
group environment. The timescale implied by red colour but visible 
spiral arms suggests that our red spirals are a likely candidate for 
the population of relative ly-faint, red g alaxies that has built up in 
clusters in recent times jPesai et alj2007h . Unless these red spirals 
are somehow 'revived', the only clear option for them is to evolve 
into galaxies resembling SOs. It is not a large step to presume that 
the objects that began their star-formation decline earliest may al- 
ready appear to have SO morphology. The existence of red spirals in 
dense environments thus lends support to the hypothesis of spirals 
evolving to SOs through environmental processes. 

It seems likely that low mass SOs do form from spirals, 
while high mass SOs were formed more directly at an earlier 
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epoch, in a manner similar to elliptica ls jPoggianti et alJ EoOll : 
iMehlert et al]|2003l : iBedregal et al]|2008l) . However, Fig. II 21 shows 
that a morphology-density relation exists at all stellar masses. 
Galaxies of all stellar masses are thus subject to environmental 
processes, though the mechanisms and the epoch at which they 
were effective may vary with mass. Candidate mechanisms, in or- 
der of declining effectiveness, include majo r mergers, minor in erg- 
ers jMihos & Hemqu ist 1994). harassment jMoore et al.l 19961) and 
cluster tidal effects ( iMerritt, 1983.) . All of these may heat or destroy 
a spiral's stellar disk, leading to an early-type morphology. Our re- 
sults clearly show that galaxy colour is transformed more readily 
than morphology, particularly for low mass galaxies, so we should 
expect to see significant numbers of red spirals, even if they will 
eventually evolve into SOs. 

Figures [TS] & [T9] provide conclusive evidence for a process 
closely related to environment which causes galaxies with iden- 
tical stellar mass and morphology to be more likely to be red in 
denser regions. This affects both spiral and early-type galaxies, 
though red spirals are more common than blue early-types. It af- 
fects low mass galaxies significantly more than those with high 
stellar mass. This process, or another, also results in an increase in 
the early-type fraction in dense environments, but a transition from 
blue to red is not generally accompanied by a transition from spiral 
to early-type. This rules out major mergers as the responsible mech- 
anism, a s they generally result in an early-type morphology (how- 
ever, seejHopkins et al. ("2008^). The most likely culprit is starva- 



tion jLarson et al. 1980i) , where the gaseous halo of a galaxy is re- 
moved by some environmental process, and therefore ceases to sup- 
ply cold gas in the disk. As the current store is consumed by star for- 
mation, there is no more gas available and star formation declines 
and eventually terminates. The main candidate mechanisms for re- 
movi ng the halo gas are ther mal evaporation jCowie & Songailal 
ram-pressure stripping jGunn&Gottlll972h . viscous strip- 
ping (Nulsen 'l982l), and interaction with the cluster tidal field 
(Merritt 1983). Direct stripping of gas from the di sks of galaxie s 
may also be important with in groups JSoseHi & Gavazzil l200i) . 
Note that at higher redshift, IPoggiantiet aL C2008) present signs 
that the environmental dependences of star formation and morphol- 
ogy may be more strongly linked than they are today. They also 
find an indication that star formation is enh anced in intermediate - 
density environments at z ^ 0.6 (se e also Bamford et al.ll200^ . 
in contrast to the independence (e.g ., Balogh et al. 2004) or sup- 
pression (e.g.. lKoopmann & Kennevll2004h of star formation with 
respect to increasing local density seen nowadays. This suggests 
that the dominant mechanism transforming galaxies in dense envi- 
ronments may change with cosmic time. 

The cessation of star formation in galaxies entering the halo of 
a larger galaxy is a key ingredient of mod els of g alaxy evolution, 
both semi-analytic models (e.g . Kauffm ann et alJI 19931) and ana - 



lytic halo models dCoorav & Sheth.2002i : [skibba & Shet"hll2008h . 
All of these models make a critical distinction between the central 
and satellite galaxies in a halo. The termination of star formation 
in satellites, and the manner in which it has been achieved, has 
generally been handled very simplistically in these models. Even 
so, with the addition of various assumptions, these models have 
met with considerable success in reproducing broad features of the 
galaxy population. Recent refinements have led to significant im- 
provements in the quantitative agreement of mode l satellite colours 
with observations jFont et ai] |2008l ; ISkibball2008h . However, there 
are many aspects which remain to be reconciled with observations. 
In particular, comparison between the morphology dependences 
predicted by these models and observations has barely begun. The 



dependence of morphology and colour on environment and stellar 
mass outlined in this paper provide clear observations with which 
to refine the next generation of models. 



5 CONCLUSIONS 

The Galaxy Zoo project has produced a catalogue of visual mor- 
phological galaxy classifications more than an order of magnitude 
larger than any previous catalogue. We have illustrated the biases 
present in this dataset, which are likely to have affected all previ- 
ous work using morphologies to some extent, but which are only 
directly quantifiable in a dataset of this size. More importantly, we 
have demonstrated a procedure for measuring and correcting for 
these biases, such that they may be reliably used for statistical stud- 
ies of morphological type fractions. Without this correction proce- 
dure the number of galaxies which may be used for such analyses 
is severely limited. 

With the de-biased Galaxy Zoo dataset we have examined 
the relationship between morphological type fractions and environ- 
ment, characterised by both local galaxy density and the distance 
to the nearest galaxy group. We reproduce the trends seen by ear- 
lier studies, but in much greater detail. The early-type fraction rises 
smoothly with increasing local density or decreasing groupocen- 
tric distance. Part of this trend is due to an increased proportion of 
massive galaxies in dense environments, which are preferentially 
early-types in any environment. However, at fixed stellar mass, a 
morphology-density relation is still present, with a similar change 
in early-type fraction with environment at all stellar masses. 

Remarkably, the morphological fraction of group members 
varies very little with the mass of the group, either determined from 
the velocity dispersion or using integrated light as a mass proxy. 

In studies based on modem surveys, galaxies on either side 
of the colour bimodality are often discussed in morphological lan- 
guage. However, colour and morphological fractions depend dif- 
ferently on environment. The fraction of red galaxies varies more 
strongly than the early-type fraction does. The result is that there is 
a considerably higher fraction of galaxies with red colour but spi- 
ral morphology in denser environments. There is also a substantial 
population of galaxies with blue colour but early-type morphology 
in low-density environments. The prevalence of red spirals peaks 
at intermediate local densities and within a few times the virial ra- 
dius of galaxy groups. These objects mostly have true spiral mor- 
phology, and are additional to the SO population. Their fraction de- 
creases in the group cores and at the highest densities, although 
this is largely due to the rapidly declining fraction of all spirals in 
these regions. The fraction of spirals that are red increases steadily 
with density. Likewise, the fraction of early-types with blue colours 
declines steadily with increasing environmental density. High stel- 
lar mass galaxies, both spiral and early-type, are significantly less 
affected by their environment than galaxies with low stellar mass. 
Most low mass galaxies, of any morphology, are blue in the field 
and red in dense environments. 

Only a small fraction of the colour-density relation is a con- 
sequence of the morphology-density relation. It is primarily driven 
by the environmental dependence of the fraction of red galaxies at 
fixed morphology: red spirals and blue early-types. Note that we 
consider morphology in terms of spiral versus early-type in this pa- 
per. It is remains possible that there is as stronger link between the 
colour-density and morphology-density relations in terms of an en- 
vironmental trend of morphology within the spiral class (i.e., Sa, 
Sb, Sc, . . . ). However, such a trend must be due to more gentle pro- 
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cesses than those which effect a full transformation from spiral to 
early-type. We therefore rule out mechanisms sufficiently violent to 
destroy galaxy disks, such as major mergers, as the primary cause 
of the colour-density relation. 

These results catalogue a complex relationship between the 
morphology, colour, stellar mass and environment of galaxies, 
which successful models of galaxy formation and evolution must 
explain. 
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Figure Al. The redshift distribution of the luminosity-limited sample anal- 
ysed in this paper (solid line) and the SDSS Main Galaxy Sample (dashed 
line) from which it is drawn. 



APPENDIX A: QUANTIFYING CLASSIFICATION BIAS 
A I Demonstrating the bias 

In the absence of evolution in the galaxy population, the true distri- 
bution of galaxy morphologies does not change with redshift. How- 
ever, the apparent brightness and size of a given galaxy depends on 
the distance at which it is viewed, and so its measured morphology 
may vary with redshift. We term this classification bias. 

For example, a spiral galaxy viewed nearby may be easy 
to classify correctly, but if viewed at higher redshift the lower 
signal-to-noise and smaller size, relative to the atmospheric seeing, 
would make the same galaxy appear as a 'fuzzy blob'. Classifiers 
with a good understanding of these issues will tend to recognise 
these limitations and assign such an object to the unclassified cate- 
gory, though not always. Non-experts, such as the vast majority of 
Galaxy Zoo participants, will more frequently classify such an ob- 
ject as an early-type. This is exacerbated by the scaling of images 
to always present objects at a similar size on-screen, as was done 
for Galaxy Zoo. However, it is likely that all morphological stud- 
ies are affected by classification bias to some extent, even including 
those studies performed by expert classifiers with the ability to vary 
the characteristics of the image display. This is the first study with 
a sample of morphologically classified objects that is sufficiently 
large for this bias to be directly quantified. 

The redshift range probed by the majority of our Galaxy Zoo 
full sample, 0.03 ^ 2 < 0.15, (see Fig. lAlt corresponds to an 
interval of 1.5 Gyr. This is 10 per cent of the age of the universe, 
during its most quiescent period. Over this period the galaxy popu- 
lation has not changed greatly. The latest studies of evolution in the 
luminosity functions of red and blue galaxies fin d a ~ 2 increase i n 
the ratio of red versus blue galaxies since z ^ 1 dFaber et alj2007h . 
which roughly corresponds to ~ 30 per cent over the redshift range 
of the Galaxy Zoo sample. However, this figure is fairly uncertain 



as it is not clear when during the z = - 1 period most of this 
evolution occurred. In demonstrating and quantifying classification 
bias in the Galaxy Zoo data we will often assume negligible evolu- 
tion. If this assumption is not valid, then the classification bias we 
measure will also include a contribution due to galaxy evolution. 
Ultimately, our aim for this paper is to correct the raw Galaxy Zoo 
classifications so that objects over the range of redshifts sampled 
may be combined. Removing any evolution, in addition to classifi- 
cation bias, will help us to achieve this aim. 

Observational limitations can introduce a variation in the pop- 
ulation of galaxies that is sampled at different redshifts. Quantities 
measured from the sample may therefore vary with redshift, even 
if the galaxy population does not. This is termed selection bias. 
For example, intrinsically faint objects can only be detected nearby, 
while luminous objects are rare and so only appear in the sample at 
larger distances where the sampled volume is larger. The apparent 
magnitude limit and volume effect dominate the variation in sam- 
ple selection with redshift. Figure [T] (in the main part of the paper) 
indicates that redshift-dependent biases are present in the GZ data, 
but does not discriminate between classification and selection bi- 
ases. However, by considering objects in narrow bins of intrinsic 
luminosity, i.e. absolute magnitude, we can remove the effects of 
this component of selection bias. 

Figure lA2l illustrates the variation in morphological type frac- 
tion as a function of redshift for galaxies in three absolute magni- 
tude intervals. If evolution is negligible, there are no selection bi- 
ases with respect to quantities other than apparent magnitude, and 
in the absence of classification bias, we would expect these type 
fractions to be constant with redshift. Note that the fractions change 
between different absolute magnitude bins due to selection effects 
and the inherent correlation between morphology and luminosity. 

Intrinsically bright objects are included in the Galaxy Zoo se- 
lection limits over a fairly wide redshift range. Their type fractions 
remain fairly constant with redshift, as expected if there is no clas- 
sification bias. However, the majority of objects with this luminos- 
ity are truly early-types, and hence free from the classification bias 
effect. The primary trend is a gradual increase in the fraction of 
galaxies classified as 'don't know', mostly at the expense of spi- 
rals. 

At more typical galaxy luminosities and fainter, a gradient 
in the type fractions versus redshift is obvious. Towards higher 
redshifts a greater fraction of objects are classified as early-type 
rather than spiral. A strong classification bias therefore appears to 
be present in the data, with apparently faint galaxies being prefer- 
entially classified as early-type, presumably because fewer details 
may be discerned given the noise. The presence of classification 
bias is also indicated by the thin lines in Fig.|2](in the main part of 
the paper), which shows type fraction trends based on the raw type 
likelihoods with selection effects removed below a given redshift 
by imposing a faint magnitude limit. 

In Fig. IA2I the tum-up in the spiral fraction at the highest 
redshifts probed in each luminosity bin is a reversal of the more 
general bias, and its origin is unclear. Plausible possibilities affect- 
ing the faintest galaxies are that the users realise they are seeing a 
blurred image and attempt to compensate with their classifications, 
that noise artifacts are more frequently interpreted as structure, or 
that early-types are being lost as a result of the effective size selec- 
tion limit. However, it does not appear to be a significant effect for 
the redshift ranges considered in the main part of this paper 

In addition to apparent magnitude, the other main factor in- 
fluencing a Galaxy Zoo classifier's ability to accurately determine 
a galaxy's morphological type is its apparent size. However, if 



© 2008 RAS, MNRAS OOO.fTltsTI 



Galaxy Zoo: morphology and colour 27 



■ early— type 
spiral 
don't know 
merger 



-23.00<M^<-22.75 
20553 galaxies 



0.1 

redshift 



o.a 



1 1 1 1 1 1 

early— t3T)e 

— — spiral 
■ ■ ■ ■ don't know 
merger 


-21.25<M,<-21.00 
45390 galaxies 




- 








0.1 

redshift 



0.2 



early-type 

— - spiral 
don't know 

— ■ - merger 



-19.75<M^<-19.50 
12316 galaxies 




0.1 

redshift 



0.2 



Figure A2. Morphological type fraction versus redshift for galaxies in three example bins of absolute magnitude from out full sample. The thick and thin lines 
con'esponds to the weighted and unweighted samples, respectively. 
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Figure A3. Ratio of the early-type to spiral measured type likelihoods as a function of absolute magnitude and physical size, for galaxies in three example 
redshift bins. The logarithmic colour scale is shown by the bar on the right. The dotted lines indicate the r = 17.77 apparent magnitude Umit, fiso^,- = 23 
mag arcsec"^ apparent surface brightness limit, and the physical scale corresponding to an angle of 1 arcsec, at the central redshift of each bin. 



we consider the measured type fractions versus redshift in bins 
of physical size (in kpc), then the apparent magnitude selection 
effects are not accounted for and dominate, togeth er with the in- 
trinsic correlation betwee n luminosity and size (e.g. lde Jondll995l : 
ICameron & Drivej2007l) . 

We therefore need to consider the measured type fractions as 
a bivariate function of both luminosity and sizj^- In this paper we 
are concerned with the early-type and spiral fractions. Studying 
the merger candidates will be an interesting topic for future stud- 
ies, but their classification is complicated by additional biases. In 
any case the numbers of objects classified as mergers are low. The 
'don't know' option provides an indication of where the classifiers 
themselves feel they are unable to determine accurate classifica- 
tions. This supplies supplementary information forjudging the un- 
biased region of parameter space, but ideally we wish to determine 
this directly from the biases in the measured type likelihoods. We 
therefore concentrate on the ratio of early-type to spiral likelihoods 
for the remainder of this paper This ratio is plotted for three ex- 
ample redshift bins in Fig. I A3 1 The number of galaxies in each 
luminosity-size-redshift bin are shown in Fig. lA4l 

Assuming evolution is negligible, and that any selection bi- 
ases are functions of only apparent magnitude and size (including 

^ We could have chosen to use surface brightness in place of either lumi- 
nosity or size, but we find that surface brightness alone is not sufficient to 
describe the behaviour of the bias. 



surface brightness, which depends on both), then in the absence 
of classification bias, we would expect the ratio of early-types to 
spirals to be constant with redshift. 

In simple terms, in Fig. |A3| changes in the colour (type ratio) 
of a particular square (luminosity-size bin) between redshift slices 
indicate a redshift dependent classification bias. This bias is due to 
the decreasing signal-to-noise and resolution for similar galaxies 
viewed at increasing redshifts. A movie is available in the online 
supplementary material, which steps through this plot in redshift 
slices and demonstrates the effect more clearly. 



A2 Constructing a baseline correction 

To quantify how classification bias changes with redshift, as a func- 
tion of luminosity and size, we first construct an estimate of the 
riei /nsp versus luminosity and size in the absence of any redshift- 
dependent classification bias. To do this, for each luminosity-size 
bin we find the lowest redshift bin that contains at least a certain 
number of galaxies, and assume that the ratio in this bin is accu- 
rate. We take the minimum number of galaxies required in a bin to 
be 30. 

From inspecting Fig. |A3| this approach appears to be reason- 
able. At low redshift the full sample includes faint, small galaxies. 
Those well inside the apparent magnitude, size and surface bright- 
ness limits should be unbiased. At higher redshifts, these limits be- 
come more restrictive, and biases grow at the faint, small end of the 



© 2008 RAS, MNRAS OOO.rTHIT] 



28 S. P. Bamford et al. 





'■. 0.0300 < 


z < 0.0400 





0.0800 < z < 0.0900 



-24 -22 -20 -IB 



-24 -22 -20 -18 




-24 -22 -20 -1 



Mr 



Figure A4. Number of early-type and spiral galaxies in each bin of absolute magnitude and physical size, for three example redshift bins. The logarithmic 
colour scale is shown on the right. Dotted lines are as in Fig. |A3l 




galaxy distribution. The scarcity of bright, large galaxies means that 
they are only seen in significant numbers at higher redshifts, where 
the survey encompasses a larger volume. However, at the redshifts 
at which n^i/nsp can be estimated for these galaxies, they are still 
well within the luminosity and size limits, helped by the intrinsic 
size-luminosity relation. On the other hand, the surface brightness 
limit may cause residual biases for the largest galaxies at a given lu- 
minosity. In order to avoid these biases, we further restrict the bins 
considered to those which are 1 mag arcsec"^ brighter than the 
surface brightness limit, 1 mag brighter than the magnitude limit, 
and with size greater than twice the angular resolution at each red- 
shift. The resulting plot should be unbiased, or as close as can be 
reasonably achieved given the available data. 

Taking this approach we obtain a baseline estimate of the un- 
biased early-type to spiral ratio versus luminosity and size at 2 ~ 0, 
shown in Fig. |A5l This figure clearly shows a region of luminosity- 
size space dominated by early-types (lower-left) and another dom- 
inated by spirals (upper-right), with a fairly sharp, curved transi- 
tion between the two. There are a small number of galaxies in the 
Galaxy Zoo full sample which are located in an area of luminosity- 
size space for which we do not have a direct estimate of the local 
early-type to spiral ratio, due to a combination of the low numbers 
of objects in these areas and their excision to avoid including poten- 
tially biased regions in the baseline estimate. In order to extrapolate 
to these regions, as well as generally removing noise and reducing 
the impact of residual bias, it is advantageous to fit a smooth func- 
tion to the binned baseline estimate. 

Motivated by the observed behaviour of the early-type to spi- 
ral ratio versus luminosity and size, after trying a variety of func- 
tions we choose the following to fit the local baseline estimate: 



risD 



Pi 



1 + exp 



f Sl(fl50 



(flso) 



+ P2, 



(Al) 



This gives a smooth step function in Mr with position and width 
varying with i?5o as 



Sl(-R50) = 

and 

S2(^?50) = n + J'2(si(i?50 



(A2) 



95 j 



respectively. 

This fit, shown in Fig. |A6| provides a smooth local baseline 
estimate against we can compare the raw early-type to spiral ratio 
versus luminosity and size as a function of redshift. In this man- 
ner we can derive a correction that will remove the majority of the 
classification bias, including all of its redshift dependence over our 
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Figure AS. Baseline estimate of the local unbiased early-type to spiral ra- 
tio versus luminosity and size. The dotted lines indicate the r = 17.77 
apparent magnitude limit, /i5o,r = 23 mag arcsec^^ apparent surface 
brightness limit, and the physical scale con'esponding to an angle of 1 arc- 
sec, at the z = 0.085 upper redshift limit of the analysis samples (see 
Sec. 12. It . The thin dashed curve indicates an apparent surface brightness of 
A'50,r = 20.25 mag arcsec"^. This approximately separates early-types 
and spirals, but clearly a simple surface brightness dependence is not suffi- 
cient to describe all the behaviour in this plot. 



sample. This method will also remove any evolution effects in this 
quantity, if present over the limited redshift range of the sample. 

To estimate the correction, C{Mr, R50, z), we simply deter- 
mine the difference between nd/risp in each redshift bin and the 
baseline just determined. This correction is plotted for three exam- 
ple redshift bins in Fig. |A7| The dependence of the correction on 
redshift is shown for a few example bins in Fig. lASI 

At this point we have a baseline correction in bins of luminos- 
ity, size and redshift, C{Mr, R50, z). A small fraction of galaxies 
lie in regions of parameter space containing insufficient objects to 
determine a direct correction. For galaxies in these bins we use the 
mean correction for the nearest neighbours within 1.2 times the 
distance to the nearest bin centre. This has the advantage of aver- 
aging over bins when there are several at similar distances. This 
extrapolation is only valid close to the region which has a directly 
determined correction, but it only needs to be, as there are very few 
galaxies far from this region. 
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Mr K 

Figure A7. Difference between the raw logiQ{n^i/ngp) in three example redshift bins and the baseline. This is the adjustment that must be applied to the 
value in each luminosity-size-redshift bin to remove the classification bias. The colour scale is shown by the bar on the right. The region enclosed by the 
dashed line is determined from the data, while outside this region the adjustment is estimated from the nearest well-determined bins. Only the region containing 
significant galaxy counts is coloured. 
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Figure A6. Fitted baseline estimate of the local unbiased early-type to spiral 
ratio versus luminosity and size. Dotted lines are as in Fig. |A5] The four 
white dots indicate the bins plotted in Fig. |A8| 



A3 De-biased Galaxy Zoo samples 

The baseline correction derived above can now be direcdy applied 
to de-bias the raw type likelihoods. To obtain the de-biased type 
likelihoods we adjust the raw likelihoods of each galaxy as 



Pcl.adj 
Psp,adj 

where 



/ V Psp / adj Pol 
1 

\ Pap / adj Pap 



(A3) 



PsP /adj 



10' 



(A4) 



and = 1 - pei - Psp- 

This effect of this adjustment is shown in Fig. lA9l The overall 
effect is to lower the early-type likelihoods, particularly for objects 
where there is some indication that the object shows spiral features 
in the raw likelihood. The largest effect is for galaxies around the 
median redshift of our full sample. For lower redshifts most galax- 
ies are well classified and need little adjustment. For higher red- 
shifts only the most luminous objects are selected, the majority of 




Figure A8. Redshift dependence of the correction to logiQ(ne;/nsp) 
for a few example luminosity-size bins, at (MrjiJso) = 
(-22.625, 7.75), (-22.375, 5.75), (-21.875, 4.25), (-20.375, 2.25), 
as indicated by the white dots in Fig. IA6I The colours of the points 
corresponds to the colour scale Fig. lA7l 



which are truly early-types and therefore do not suffer from the 
bias. 

The success of our de-biasing procedure is demonstrated by 
Figs.[2]&|7]in the main body of this paper, which show that the 
method works as expected. 



APPENDIX B: 'FINGERS OF GOD' CORRECTION 

The high velocity dispersions in galaxy groups lead one to grossly 
overestimate the line-of-sight distances between the group centre 
and its member galaxies, when a cosmological conversion between 
redshift and distance is assumed. We correct for this effect by re- 
ducing the distance to the group centre for galaxies with redshifts 
consistent with the group redshift and velocity dispersion, and close 
to the group in projected distance. Specifically, we determine the 
line-of-sight comoving distance between the galaxy and group cen- 
tre as 
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Figure A9. De-biased type likelihoods versus the raw values direct from 
the Galaxy Zoo classifications. The dotted diagonal line indicates the one- 
to-one relation. Points above and below this line refer to the spiral and early- 
type likelihoods, respectively. The solid and dashed lines trace the median 
and quartiles or the points. 



dins — dc 



= 1/0^ 



-C{max(i?v 



,(B1) 



where A, B, C and D are chosen to ensure the following be- 
haviour. For galaxies with projected distances ^ iivir: dios = 
when 2gai = Zgrp; dios = RvIt: whcn l^gai — ZgrpI — C,', and the 
adjustment is less than 1 per cent when [zgai — Zgrpl ^ 2(^. Here, 
^ = So-grp (1 + Zgrp)/co is the redshift difference corresponding to 
three times the group velocity dispersion, dc is the comoving dis- 
tance between the galaxy and group redshifts assuming only cos- 
mological motions, and rfpro is the projected distance between the 
galaxy and the group's brightest member assuming they both lie 
at the group cosmological redshift. For galaxies with dpro > -Rvir 
the adjustment is decreased to 90 per cent of its ^ i?vir value by 
dpro = 2i?vir and 1 per cent of its ^ R^ir value by dpro = 5i?vir. 

We furthermore normalise the corrected distances by 7?vir of 
the nearest group to account for scaling of any potential influence 
with group size. Figure lBTI illustrates the effect of the correction for 
model galaxies in a typical C4 group. A comparison of the redshift- 
space distances, dc4 ~ ydl + df^ and the normalised, corrected 

distances, Dc4 = \/ df^s + '^pro / ^vir, is shown in Fig. lB2l 
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Figure Bl. Illustration of the 'Fingers of God' suppression method. 
Coloured lines show the relation between redshift and the logarithmic dis- 
tance, in virial radii, to a group with properties typical of those in the C4 
catalogue: z = 0.06, cr = 680 km s^^ and iJvir = 0.74 Mpc. Thin, solid 
lines plot the relation assuming no peculiar velocities. In this case, small 
departures from the cluster redshift imply large distances. However, in re- 
ality, peculiar velocities impart significant redshift differences on galaxies 
that are actually members of the group, and so at small distances from one 
another. Thick solid lines show the relation with the 'Fingers of God' sup- 
pression method, which attempts to account for this. The red, orange, green 
and blue lines indicate the effect on galaxies at projected separations of 0. 1 , 
0.5, 2, and 3 -f?vir. respectively. 
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Figure B2. Redshift-space distance of each galaxy to its nearest C4 group, 
dc4. versus a measure of tliis distance incorporating a correction for the 
'Fingers of God' effect and normalised by virial radius, Dc4- Contours 
enclose 50, 75, 90, 95 and 99 per cent of galaxies in the luminosity-limited 
sample. 
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